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Abstract 

The quantum double of the Haagerup subfactor, the first irreducible finite 
depth subfactor with index above 4, is the most obvious candidate for exotic 
modular data. We show that its modular data PHg fits into a family T>'^ilg2n+i, 
where n > and cj E Z,2n+i- We show P''Hg2„+i is related to the subfactors 
Izumi hypothetically associates to the cyclic groups Z2n+i- Their modular data 
comes equipped with canonical and dual canonical modular invariants; we com- 
pute the corresponding alpha-inductions etc. In addition, we show there are 
(respectively) 1, 2, subfactors of Izumi type Z7, Zg and Z|, and find numer- 
ical evidence for 2, 1, 1, 1, 2 subfactors of Izumi type Zn, Z13, Z15, Z17, Z19 
(previously, Izumi had shown uniqueness for Z3 and Z5), and we identify their 
modular data. We explain how PHg (more generally T>^iig2n+i) is a graft of 
the quantum double T>Sym(3) (resp. the twisted double T>^ D2n-\-i) by affine 
so(13) (resp. so(4n^ ->rAn + 5)) at level 2. We discuss the vertex operator alge- 
bra (or conformal field theory) realisation of the modular data P'^Hg2„_|_i. For 
example we show there are exactly 2 possible character vectors (giving graded 
dimensions of all modules) for the Haagerup VOA at central charge c = 8. It 
seems unlikely that any of this twisted Haagerup-Izumi modular data can be 
regarded as exotic, in any reasonable sense. 
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1 Introduction 

From the early years of conformal field theory (CFT) (see e.g. [511 [60]), we find: 

Speculation. The standard constructions (orbifolds, cosets, simple- current exten- 
sions, ...) applied to the basic theories (lattice compactifications, affine Kac-Moody 
algebras, ...) exhaust all rational theories. 

Perhaps this should be applied to the modular tensor category rather than the 
full CFT, in which case this would constitute a sort of generalised Tannaka-Krein 
duality holding for modular tensor categories, where the dual to the category (i.e. 
the analogue of the compact group) is (say) a vertex operator algebra constructed in 
a standard way from the basic examples. 
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Possible counterexamples have been around just as long. For instance Walton [59] 

explained that rank-level duality applied to conformal embeddings of affine algebra 
CFTs (the simplest being Ai iq C 6*2,1, yielding an extension of ^9,2) yield what seem 
to be new CFTs. The possibility that rank-level duality is somehow inherently sick 
was eliminated by Xu [61], who realised many of these examples with completely 
rational nets of subfactors (these should correspond to rational CFT). Perhaps all 
that this accomplishes is to insist that rank-level duality (or in subfactor language the 
mirror) should be included as one of those standard constructions. More recently, [12] 
proposed four modular data (each with 6 or 7 primaries) as possible counterexamples 
to this speculation, but didn't show they could be realised by a rational CFT (or 
vertex operator algebra or net of subfactors). 

Perhaps the most obvious place to look for a truly exotic example is the Haagerup 
subfactor [T]. It is the first irreducible finite depth subfactor with index greater than 
4 — it has index (5 + a/TS) /2 ^ 4.30278. It is generally regarded as exotic, since (so 
far) it can only be constructed by hand without any natural algebraic symmetries. 
However it is nonbraided and so to get a braided system (which has a chance to 
correspond to the fusion ring and modular data of a rational CFT) we should take its 
quantum double, or equivalently asymptotic inclusion (see Ocneanu as in [251 Chapter 
12]) or Longo-Rehren inclusion [l9l|50]. Its modular data 2^Hg was computed in [38] . 
and subfactor realisations of its modular invariants were studied in [26] . 

Question 1. Is this Haagerup quantum double realized by a rational CFT (or rational 
vertex operator algebra — see Definition 3 below — or completely rational net of 
subfactors l4^), in the sense that they share the same modular data T>B.g? 

We would want a deeper relation between the Haagerup double and the corre- 
sponding vertex operator algebra (VOA) than merely that their modular data coin- 
cide, but for now that can suffice. 

A more direct realisation of the Haagerup subfactor C M in a rational CFT 
could be that its A^-A^ sectors (say) form the algebra of defect lines (or full system) 
of a rational CFT. After all this will in general fail to be commutative. However, 
even this cannot happen: the modular data of the quantum double of a full system 
must be in factorised form X ® This is because the quantum double of the 

full system is equivalent to the double of the fusion algebra X by [52], [531 Cor. 2.2], 



and this double factors as X ® X°^^ by [24^ Prop. 2.2] as long as the braiding on 




X is non-degenerate. However, the modular data of the double of the Haagerup 
is easily seen to not be in factorised form [35]. There remains the possibility, which 
we won't explore in this paper, that the Haagerup systems are the full system of a 
degenerately braided system. Note that much of the theory of modular invariants, 
including alpha-induction and the full system, holds for degenerately braided systems 



Question 2. IfVHg is realised by a rational VOA, say, then is that VOA exotic, in 
the sense that it cannot be constructed from standard methods and examples? 

The three known constructions of the Haagerup subfactor (namely Haagerup's 




[mini El El HO]. 
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connection computation [3ll[T], Izumi's Cuntz algebra construction [38] and Peters' 
planar algebra construction [51]) are largely combinatorial tours de force (although 
Izumi's suggests an underlying cychc group), and for this reason the Haagerup subfac- 
tor is generally regarded as exotic. Even given this, it is conceivable that its quantum 
double "PHg could be constructed directly in more standard ways. [35] confirm that 
PHg does not fall into the simplest possibilities, namely the modular data of an affine 
algebra, lattice, or finite group. 

We argue that the answer to Question 1 is yes, and that for Question 2 is no. 

In particular, we explain why such a VOA should be a conformal subalgebra of 
the central charge c = 8 VOA V{EqA2) corresponding to (root) lattice Eq © A2 (by 
conformal subalgebra we mean a sub VOA of identical central charge and identical 
conformal vector). A more familiar conformal sub VOA of V{EqA2) is an order-2 
orbifold, realising the modular data of the quantum double of the symmetric 
group Sym{3) = S3. An early indication for us that the Haagerup is related to S3 
was the similarity of their tube algebras. Tube algebras were introduced by Ocneanu 
(see e.g. [25]) as a way of computing the irreducible objects and fusion rules of the 
quantum doubles of subfactors, and developed further by Izumi [371 ESj in particular 
in analysing and determining their modular data. The list of modular invariants and 
nimreps for the doubles of the Haagerup and 5*3 [T51 EZ] are also strikingly similar, 
and in fact the modular invariants lead us to Eq © A2. The connection between 
modular data and VOAs is made much more explicit using character vectors. A 
second, equally promising possibility for a VOA realisation of I^Hg, by GKO cosets 
using the affine algebra VOA V(i?6,2), is discussed but not analysed in detail. 

We generalize the modular data I^Hg in two directions, by fiting it into an infinite 
sequence, and showing the nth term in this sequence can be twisted by Z2„+i. The 
role of 5*3 is now played by the twisted quantum double of D2n+i, and V{EqA2) 
(which realises the VZ3 modular data) by a holomorphic orbifold by Z2n+i. Subfactor 
realizations of V^}lg2n+i should be provided by Izumi's hypothetical generalisation of 
the Haagerup: we construct a unique subfactor of Izumi type which realises V^Hg^ for 
all 1 < z/ < 19 and conjecture that indeed this continues for every n. For this reason 
we suggest calling D'^Hg2„_,_i Haagerup-Izumi modular data, von Neumann algebra 
realisations of twisted Haagerup-Izumi modular data P'^Hg2„+i are still unclear. 

Question 3. Can we orbifold a VOA by something more general than a group? 

If so, this should provide a simple construction of VOAs realising P'^Hg2„_,_]^, 
starting from a holomorphic orbifold by Z2n+i- Perhaps this is not that unreasonable 
— after all, a subfactor itself is really a generalised orbifold, in the sense that the 
quantum double of a subfactor M*^ C M of fixed points recovers VG. Indeed the 
whole motivation of subfactor theory (see e.g. [25]) is to understand non-grouplike 
quantum symmetries through subfactors or inclusions of factors C M, and asso- 
ciated notions of N-M, N-M, M-M irreducible bimodules or sectors, their fusion 
rules, paragroups, A-lattices or planar algebras. This notion of generalised orbifold 
construction also features in the framework of |30j . 

Our work suggests a further generalisation. Consider a quadruple {K, H, a, u) 
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where K and H are finite groups, and H acts freely on i.e. tfie group homo- 
morpliism a : H ^ Aut(-ft') obeys a{h)k = k iS h = e or k = e. Freeness implies 
tlie projection h i— Out(-ft') is injective; such a semi-direct product D = KxiaH is 
called a Frohenius group. The twist uj G H^{Bk', T) should be compatible with a in 
the sense that it lies in the image of the natural map H^{Bd', T) — )■ H^{Bk', T). The 
Haagerup subfactor is associated to (Z3,Z2,a(l) = — , [0]) and D = S^. 

Wildly Optimistic Guess. Let {K,H,a,uj) be any quadruple defined above, cor- 
responding to Frobenius group D = Ky\aH. There is an irreducible finite-depth 
subfactor associated to the pair {K, H); the triple {K, H, a) is realised by a Q-system 
of endomorphisms. To any such quadruple there is a rational VOA V{K,H,a,u), 
and a completely rational net of subf actors realising the modular data of this twisted 
quantum double. This VOA is a generalised orbifold (controlled in some sense by 
H) of a holomorphic orbifold Vk by K. This VOA Vk o,lso contains a holomorphic 
orbifold Vd; both V{K, H,a,u) and Vd contain a common rational VOA. For K an 
odd cyclic group and H = acting by n ^ —n, this recovers the modular data given 
in Section 3.2 and also a subfactor of Izumi type K, i.e. a solution to equations 
(7.1)- (7. 5) of lEBI; the K = special case recovers the Haagerup subfactor. 

This guess suggests that the Haagerup subfactor itself cannot be regarded as 
exotic. But this guess is in some control only for the odd cyclic case discussed above, 
where it seems quite plausible. The reason for requiring freeness is that this condition 
plays an important role in the derivation of V^Hg^^ in Section 3.2, as well as in [38] : 
probably it can be dropped but at a cost of complexity. 

The possibility that the Haagerup subfactor is related in some way to Ss appears 
to have been first made in [21]. The present paper clarifies, deepens, and generalises 
that relation. The remainder of this first section recalls modular data and modular 
invariants (Section 1.1), reviews the basic theory of subfactors and the tube algebra 
(Section 1.2), and discusses the modular data and character vectors of VOAs (Section 
1.4). We also explain in Section 1.3 how to recover the original subfactor from its 
quantum double — unlike the rest of Section 1, this subsection contains original 
material. Section 2 develops the relation between the Haagerup subfactor and 5*3. 
Section 3 puts the Haagerup modular data into a sequence and twists it. We also 
introduce the notion of grafting and discuss the role of affine so(z/^ + 4) at level 2. 
Section 4 relates our untwisted sequence to Izumi's hypothetical family of 3-star, 5- 
star, 7-star, ... subfactors. Section 5 discusses VOA interpretations for our twisted 
sequence, and explains the Haagerup-dihedral diamond which generalises Section 2. 

The Haagerup subfactor arose in Haagerup's classification [5^ of irreducible finite 
depth subfactors of index between 4 and 3 + y/S ~ 4.73205. Two other subfactors 
N (Z M (both also regarded as exotic) appear there: the Asaeda-Haagerup subfactor 
with Jones index (5 + a/I7)/2 ^ 4.56155 [1], and the extended Haagerup subfactor [1] 
with index ^ 4.37720. Both seem unrelated to our family, and neither lies in a known 
sequence. As with the Haagerup subfactor, neither their A^-A^ nor M-M systems are 
braided. Of course Questions 1 and 2 should also be asked of them, but no serious 
work can begin until the modular data of their doubles has been computed. 
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1.1 Modular data and modular invariants 



Modular data arises naturally in several contexts (for instance CFT, subfactors, and 
VOAs) — see [HT] for a review. In CFT, modular invariants correspond to the 1-loop 
closed string partition function. 

Definition 1. Modular data ($, 0, S, T) consists of a finite set $, an element G 
and matrices S = (S'jj)jjg$ and T = (Tij)ij^<;,, such that: 

a) S, T are unitary, S is symmetric, T is diagonal; 

b) S"^ = (ST)^, = I , the identity matrix; 

c) C := S"^ is a permutation matrix, and row of S consists of nonzero real numbers; 

d) for all i,j, /c G $, the following quantities are nonnegative integers: 

We call modular data ($, 0, S*, T) and ($', 0', 5*', T') equivalent if there is a bijection 
: $ ^ $' such that 0(0) = 0', S'^^^^^^^^^ = Sij and T;(,),<^(,) = T,,, for all i,j G $. 
The i G $ are called primaries and G $ is called the vacuum. The order 1 or 
2 permutation C is called charge- conjugation. Tq q = e~'^ic/i2 ^ number c 
(only defined mod 24 by the modular data) called the central charge. Equation 
(11. ip is called the Verlinde formula; the quantities Nf^ are called fusion coefficients 
and comprise the structure constants of a commutative associative algebra called the 
fusion ring. The coefficients iVi^-fc := Nf^J are invariant under all 6 reorderings of 
indices. Also, N^^ij = Si^cj- 

Condition b) says that modular data generates a representation p of the modular 
group SL2(Z) through the assignment 

s=p(° .r=p(j ;). (1.2) 

Most examples considered in this paper have the additional property that C = I, i.e. 
S is real and the representation is factors to PSL2(Z). We assume throughout that 
the first row of S is strictly positive — such modular data is called unitary. 
From Definition 1 we obtain 

Sij = Ti^iTjjTQ Q Tk^kSkfiN^j . (1.3) 
fee* 

Thus if T is known exactly but S only approximately, (ll.3p together with the integral- 
ity of the fusion coefficients (11. ip can be used to determine 5* exactly (the quantum- 
dimension Si^o/Sofi is the Perron-Frobenius eigenvalue of the matrix Ni = {Nf'j), and 
5*0,0 > is determined by the square ^j('S'j,o/'S'o,o)^ of the global dimension I/Sqa). 
Equivalently, if S,T and S',T are both modular data sharing the same T, and S' is 
sufficiently close to S, then S = S'. This is how we'll identify the modular data of 
the double of subfactors in Section 4.1, given numerical estimates of S. 
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One consequence of tliis definition is that the entries of S and T he in a cy- 
clotomic extension Q[C,l] of the rationals (throughout this paper we write for 
exp[27ri/A;]). For any Galois automorphism a G Gal(Q[^L]/Q) ^ Z^, ($, 0, ct5, ctT) 
wiU be (generally inequivalent) modular data with identical fusion coefficients, where 
a acts entry- wise on S and T. It can be shown [Hj that for any a G Gal(Q[^i;,]/Q), 
there is a permutation j i— )■ aj of $ and choice So-(j) of signs such that 

a{Sij) = s^{i) S^ij = s^{j) Si^„j . (1.4) 

An important class of examples of modular data comes from finite groups G. 
Write e for the identity and g'^ for h~^gh. Select a cocycle u G Z^{Bg]T), where 
we write T for the unit circle in C; the corresponding modular data depends up to 
equivalence only on its class in H^{Bg', T). Define 

e,{h, I) = uih, g\ I) uj{g, h, I) uJ{h, /, • (1-5) 

Then the restriction of each 6a to the centraliser Za := {h E G : ha = ah} in G is 
a normalised 2-cocycle. In this paper we are only interested in the cohomologically 
trivial case where all 6a are coboundaries (this will happen for instance whenever the 
Schur multiplier H^{Bz^', T) of each centraliser Za is trivial). This means that there 
exist 1-cochains : 2^ -> T for which ea{e) = 1 and both 

Saih, g) = {5ea)ih, g) = ea{h) ea{g) ea{,hg) , (1.6) 
(^x-^ax{x~^hx) = 6a{x, x'^hx) 6a{h, x) ea{h) , (1.7) 

for all g,h E Za and x E G. The primaries are all pairs (a, vr) where a runs through 
representatives of conjugacy classes in G and n E Za = Irr(Za) are irreducible repre- 
sentations {irreps). The vacuum is (e, 1). The modular matrices are 



Sia,x),(a',x') = ly My I xig^' 9 ^) X (o,^) ^aiga' g ^) ega' g~i (a) , (1.8a 

' 3eG(a,a') 

x(q 

X(e 



^{a,x),(a',x') - ^«X^x,x'T-77Y ^«(") ' (1.8b) 



where G{a^a') = {g E G : a^b = ba^}. In the following we let V^G denote this 
modular data. See e.g. [121 [IZJ [13] for more details; we recover this formula next 
subsection using tube algebras. 

Another class of modular data is associated to even positive-definite lattices. Let 
L be such a lattice, and n its dimension. Then $ = L*/L where L* is the dual of L. 
The vacuum corresponds to the coset [0]. Then 

S[uUv] = \L*/L\-^/^exp[27nu ■ v] , (1.9a) 
T[u],lu] = exp[7riM ■ u — 7iin/12] . (1.9b) 

The final class we need is associated to affine nontwisted Lie algebras g*-^-* and 
positive integers k (see e.g. |12])- Let g^k denote the corresponding modular data: 
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$ is the set P^{q^^^) of integrable highest- weights of level k, the vacuum is /cAq, and 
explicit formulas for the matrices S and T are given for instance in Ch.l3] or 
PT| . Finite group, lattice and affine algebra modular data (and as we'll see next 
subsection, the double of a subfactor) is always unitary. 

Definition 2. A matrix Z = {Zij)ij^<i, is called a modular invariant provided 

i) all entries Zij are nonnegative integers and Zq^ = 1; 

ii) ZS = SZ and ZT = TZ . 

Examples are Z = I and Z = C. A modular invariant is often written in equiva- 
lent form as the formal expression 

Z = Zijchichj . (1-10) 

Unitary modular data has only finitely many modular invariants. The Galois sym- 
metry ( !Li|) implies, for any a G Gal(Q[^L]/Q), 

= s,(2)s,(j)M,,,,,, (1.11) 
^ implies s^{{} = s^U) . (1.12) 

1.2 Subfactors, quantum doubles and tube algebras 

We refer to [HI for the basic theory of subfactors, principal graphs etc, [IS] for 
the theory of sectors, and [HI O HO] for the theory of alpha-induction. 

Given a type III factor N, let N'^n denote a finite system of endomorphisms on 
N [m Defn 2.1]. Write I1{n'^n) for the endomorphisms which decompose into a 
finite number of irreducibles from N'^n- For A,p G T,{nXn), the intertwiner space 
Hom(A, p) is a finite-dimensional Hilbert space. Write (A, p) for its dimension. The 
sector [A] identifies all endomorphisms Ad(M)A for unitaries u in the target algebra. 

Suppose N'^N is nondegenerately braided [HI Section 2.2]. Among other things 
this means A,/i G n^n commute up to a unitary e = e(A,/i), i.e. Ap = Ad(e)/iA, 
and these unitaries {e(A,yu)} can be chosen to satisfy the braiding-fusion relations. 
Unitary modular data S,T is obtained from this set-up by the intertwiners asso- 
ciated to the Hopf link and twist [33 EH]; here $ = n'^n and 1 is the identity 
endomorphism, and the corresponding fusion ring is realised by composition in N'^n'- 
[A][/i] = E<m[^] where A^f^i^^j = ([A][/i], [u]). 

Modular invariants are recovered as follows. Suppose we have a subfactor A^ C M. 
Let L : N ^ M he the inclusion and I : M N its conjugate. Then 9 = II is called 
the canonical endomorphism and 7 = a its dual canonical endomorphism. Suppose 
9 is in S(ArAfAr). Using the braiding e"*" := e or its opposite e~ := e"^, we can lift 
an endomorphism A G N'^n of A^ to one of M: := 7~^Ad(e^(A, 6'))A7. Then 
:= {a'l,a~) is a modular invariant [HI |20]. The induced a^(Ar^Ar) = M'^m 
generate the full system m'^m- The N-M system n'^m (resp. the M-N system 
M'^n) consists of all irreducibles in l9"' (resp. I7"). By the nimrep we mean the nXjsi 
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action on the N-M system N'^m- This is one of 8 (6 independent) natural products 
pXq X qXr — )■ pXji among the sectors, one for each triple P,Q,R& {M, A^}. 

Many examples (e.g. the Haagerup and finite groups) arise naturally as nonbraided 
systems of endomorphisms. To get a braided system, one takes the quantum double 
(asymptotic inclusion) of N'^n- This can be realised by the Longo-Rehren inclusion 
A G B where B = N ^ iV°pp (see e.g. |l9l |37]): we are interested in the double 
©(^tAtv) system on A (which contains the A-A system of this subfactor, cf. Remark 
(i) p. 146 of [2Z|) whereas the full system b'^b is simply j^X^^ ^X^"'^'^ (see Theorem 
1 below) with dual canonical endomorphism = T^^^j^Xn^, ® 

If is an M-M system of a subfactor C M, then the dual Longo-Rehren 

inclusion is "T: Ai C M ® M°pp =: Bi, where 7 = 'Lt= XlrjeM^'M ^ ® r]"^^. Here Bi'^Bi 
is M'^M ® M'^M^ but we can and will identify the double V{]^Xn) on A with the 
double T>{mXm) on Ai, and Ai with A. 

The induction-restriction graph of the quantum double system is constructed in 
the following way |3^. First, given p G T,{^X^), a system of unitaries {£p{C)} 
is called a half-braiding of p if Sp{C,) G Hom(p,^,,^p) for some ^ G n'^n, and 
xSp{() = ^{£p{vi))£p{^)p{x) for every X G Hom(C,,^?7), as in [371 Defn. 4.2]. This 
can be described by a matrix representation for an orthonormal basis {wp(0«} of 
Hom(^,p), if we set 

^p'(0(..),(C,.) ■■= Pdwp{0*)£l{i)wp{^). G Hom(7^e,eC) • 

Then the even vertices of the quantum double system are labelled by inequivalent 
half-braidings and the odd vertices are labelled by (^ ® id"''^)/, with ^ running in 
nX]si. Finally vertices £p and (^ ®id°''P)t are joined by (p, ^) edges, see [STJ page 154]. 
The principal graph of A C -B is the connected component containing (id ® id°''^)<.. 

Incidentally, the forgetful functor sending the half-braiding to the representa- 
tion p is an algebra homomorphism, and indeed has a name: it is alpha-induction. 
See also the discussion at the end of the next subsection. 

The tube algebra is an effective way to compute the modular data 5, T of the 
quantum double. In Section 2 we compare the tube algebras of the Haagerup sub- 
factor and the group 5*3. Their similarity is the first indication of a possible relation 
between them. The tube algebra of n^n is the finite dimensional C*-algebra 

Tube(^A'^)= Hom(e-C,C-^)- (1-13) 

The simple summands of Tube(jvAf7v) are labelled by inequivalent half-braidings £^. 
The modular data S", T is obtained from the half-braidings by (2. 8), (2. 9) of [38] . 

Let's consider in more detail the special case of tube algebras of finite groups. Let 
G be a finite group with identity e. Given a type HI factor A^, write Int(iV) for the 
group of inner automorphisms and Out(A^) := Aut(A")/Int(iV). Let a : G — )■ Out(A^) 
be a homomorphism. We take a lift a : G — >■ Aut(A^), as in [57] . 

Since [ag][a/i] = [oigh] as sectors, there exists a unitary Ug^h € N satisfying 
Oig-oih = A.d{ug^h)oigh- In particular Ug^h £ Hom(Q!g/j, ag-at). Associativity {agah)(yk = 



9 



ag{ahak) implies Ad{ug^h)ugh,k ■ aghk = M{ag{uh,k)ug^hk) ■ ctghk- Then since is a 
factor, there exists a scalar u{g, h,k) satisfying Ug^hUgh,k = ^{9^ h, k)ag{uh,k)ug^hk^ 
i.e. u is an element in Z^{Bg'i"^)- Conversely, every element in Z^{Bg'i"^) arises as 
such an u [FTlliD] . 

Set c{g,h) := UhghU* i^. The product and *-structure of the tube algebra 
Tube(ArrYAr) are: 

c{g, h) c{k, I) = 5k,gH Ogih, I) c{g, hi) , (1.14) 
c{g, h)* = co{gh, h-\ h) u{h-\ gh, h-') uj{h-\ h, c{g\ h'^) (1.15) 

(recall (11. 5p ). Consequently, 

c(e,fif)c(e, /i) = c(e,5(/i) and c(e, ^f)* = c(e, gf"^) , 

and thus the group algebra C[G] of G is a C*-subalgebra of Tube(ArA'Ar). The identity 
of the tube algebra is ^g^Qc{g, e). 

Let Kg = {g'^ : h G G} denote the conjugacy class of an element g ^ G. Then 
we have the decomposition Tube(7vA7v) = ©x„ Tube(ii'„) where the sum ranges over 
the conjugacy classes of G, and Tube(i^^,i) = 0gg^ i^QB.om{agai,aiagi). Hence 
C[G] = Tube(i^e)- 

This decomposition can be further refined. Consider first a trivial twist u 
The outer action a gives a subfactor N G N xiG = M. The A^-A^ system is {a^} = G 
whereas the M-M system is the irreps G. The tube algebras Tube(G) ^ Tnhe{G) are 
Morita equivalent as finite-dimensional C*-algebras [H7] , so their centres, the quantum 
doubles of G and G, are identified. The simple components of Tube(G'), equivalently 
the even vertices of the quantum double of G, are labelled by pairs {Ku, vr) where 
TT G Zu are irreps of the centralizers Zu, see e.g. or [221 Sect. 4]. To every such 
pair {Ku,7t) is attached the endomorphism P{Ku,n) £ ^{G) defined as 

P{Ku,7r) = dim(7r) ^ at ■ (1.16) 

In the dual system G, the endomorphism P(Ku,tt) ^(^) attached to such a pair 
{Ku, vr) is the Mackey induction Ind^^vr from Zu to G, decomposed into irreps of G. 

The number of inequivalent half-braidings associated to P{Ku,-k) equals the number 
of inequivalent irreps vr' G Zu such that dim(7r') = dim(7r). The half-braiding for 
P = P{Ku,n) in its matrix decomposition is 

£;i0^u.,uu^,) = Mpr\-')id^, (1.17) 

for all G p^^ZuQ, u^, G A'„, where i,j label the basis vectors in the representation 
space of TT. 

All this extends to the twisted case. The 6a s of (II. 5p are normalised twisted 
cocycles on G, namely from (I1.14p they satisfy 

Oaix, y) Oaixy, z) = 9aix, yz) 9^-iax{y, z), yx,y,z eG. 
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They thus describe projective representations of Za- As we know [T7], the primary 
fields in the model twisted by a given 3-cocycle uj consist of all pairs (Kai ^) where 
now TT e 9a-lrr{Za). 

In the cohomologically trivial case discussed in Section 1.1, this is immediate: we 
can twist the formula f ll.l7p for half-braidings by e; inserting into (2. 8), (2. 9) of [38] 
recovers (11.81) . 

1.3 The canonical and dual canonical modular invariants 

This subsection contains original material. 

Suppose we have a subfactor C M, with A^-A^ system denoted by A and M-M 
system A. The double VA = LR{A) of A is realised as sectors on a factor A in 
A^ (8) A^°PP =: B. Alpha-induction reverses this. That is, if you take the appropriate 
modular invariant in the double, and take alpha-induction from A to B, then the 
full system is A (g) A°pp. In some sense then, the factor A and the modular invariant 
Za remembers the original system A on A. In fact it is a consequence of Theorem 1 
below that the A-B system is just A regarded as a nimrep. 

To get B, only the canonical endomorphism 6* on A is needed (because of Theorem 
1 below, we can read 6 off from the modular invariant), plus the associated Q-system 
on 6. That is, all the information about B and the full system A ® A°pp is carried 
in the primary fields and modular data of the double, the modular invariant Za and 
its vacuum block 6. We suggest this modular invariant Za be called the canonical 
modular invariant. 

Dually, there is an inclusion A G M ^ M°pp =: B, where the dual A acts on 
M, but the same double acts on A. This comes with another modular invariant Z^, 
with a corresponding 6 which regains A ® A°pp as the full system. Here (again from 
Theorem 1 below) the A-B system is just A. Again the modular invariant Z^ — the 
dual canonical modular invariant — and 6 together encode the original system A. 

This means we should regard the double VA = T>A as equipped with two canon- 
ical modular invariants Za and Z^, from which using say alpha- induction one can 
recover the full system B-B and the nimrep A-B entirely in terms of A or A. In 
particular the full system is A® A°pp or A® A°pp, depending on the choice of modular 
invariant Za or Z^, and the nimrep is A or A. 

Note that in both cases the A-B system here is an algebra. For most modular 
invariants the A-B system is only a module, but when the inclusion is type I, which 
implies for instance that the modular invariant is a sum of squares, the A-B sys- 
tem will necessarily be an algebra (A C -B is type I iff the equivalent conditions of 
Proposition 3.2 in [5] hold; in the nets of subfactors setting A(J) C B{I), this means 
the extended net is local). This suggests that both inclusions A G N ® ]\f°pp and 
A G M ^ M°PP are type I. In fact more is true: 

Theorem 1. The inclusions A G N A^°pp and A G M ® M°pp are both type I. 
Moreover, we have {ZA)i,j = (ZA)i,o(-^A)i,o ; ^-e. Za = \ J2i^i,o'^^if (similarly for 
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Z^). The vectors u and v with entries Ui = (^A)j,o one? Vi = (Z^)j^o ore eigenvectors 
with eigenvalue 1 of both S and T . 



Proof. Corollaries 6.3 and 6.4 of [3Z] show that, at the level of subfactors, B'^b ■~ 
a+iB^B) C A (g) 1 and b^^ := o'Ia'^a) ^ 1 ® A°pp. The neutral system b^b = 
B^B ^ B'^B is ^^^s contained in A ® 1 fl 1 ® A°pp and hence equals 1. But dim B'^b 
is computed in [H Prop.3.1], for such a nondegenerate system like ours, to be 

>J0,0 — >-'o,o — >-'0,0 

which matches the dimensions of A (g) 1 and 1 ® A°pp. Therefore we get the equalities 
bX^ = a (g) 1 and = 1 ® A°pp. 

Any canonical endomorphism 9 is bounded below by the vacuum sector: 9 > 
J2\^>^,o^- This inequality is a special case of {9\,ii) > (a^l',^^) ([121 eq.(37)]; see 
also IP Thm.3.9]). But dim 6* is given by 

dim 6' = dim(a) = dim(il) = ^ dim(^ ® C^^) = ^(dim^^ 



'S'0,0 

agreeing with fll.lSp . and thus 9 = XIa ^^^o-^- By [81 Prop. 3. 2], this implies the 
inclusion A C -B is type I. 

Because of alpha-induction, we know a modular invariant takes the form Z = 
X x^r,fi, where r runs over all sectors of the neutral system and bf^^ = (r, a^) 
are the branching coefficients. Because this system is type I, we have 6"^ = b~ . 
Because the neutral system is trivial, we have only one r. Thus the modular invariant 
Za corresponding to 9, and Z^ corresponding to 9, take the desired forms. The 
statements about the vectors u, v follow from modular invariance. QED 

From the theory of alpha-induction, we know that A (respectively A) is non- 
commutative iff some {Z/s)i j > 1 (resp. some {Z^)i j > 1) — see [HI Cor.6.9].[T2| 
Thm.4.11]. For type I inclusions, we know from Section 4.1 of [7] that the A-B system 
A^B is isomorphic to b^b- 

We will call any modular invariant of the type Z = \ Zi^ chi\'^ a monomial 
modular invariant. 

Of course all this applies for a system A of endomorphisms not necessarily coming 
from a subfactor — we expect this is relevant for our twisted Haagerup data. Again 
A is recovered from the double PA and the canonical modular invariant Zj\, which 
again is monomial. The only difference is that there is no A. 

Suppose a finite group G acts by outer automorphism on a type III factor A^. Then 
as mentioned last subsection, by taking a crossed product we have a subfactor N G M 
where the A^-A^ and M-M systems are identified with G and G respectively. Taking 
the Longo-Rehren inclusion A C A^® A^°pp, the doubled A- A system can be identified 
with the untwisted quantum double VG. We will sometimes find it useful to use the 
K-theory language developed in [221 [23] . which identifies T>G with the equivariant 
K-group Kq{G) = K\_^{G x G) where G acts on G by the conjugate action and 
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A = A{G) = G acts diagonally on the left and right of G x G. The neutral system 
is trivial {K^{e, e) = Z), with sigma-restriction aid = ^vg — J2xeG dimx (e, x) ■ The 
full system is A x A°pp = K^{G x G) with canonical modular invariant 

2 



e) ch 



(e,x) 



:i.i9) 



the sum over x G G. The dual LR- inclusion is A C M®M°pp, where the doubled A- A 
system is again the quantum double of G, the full system is G x G°pp = Kq^q{G/G x 
G/G) and the neutral system K^{G/G x G/G) is trivial. Here sigma-restriction is 
given by aid = 9dg — ^(S'jid) where g runs over representatives of all conjugacy 
classes, and the dual canonical modular invariant is 

^A= |E^V,i)f ' (1-20) 

where the sum runs over representatives a of all conjugacy classes of G. We re- 
cover from fll.19p . fll.20l) that A = G is always commutative, whereas A = CG is 
commutative iff all dimensions = 1, i.e. iff G is abelian. 

By the Galois theory of [3Z], a subgroup K < G induces an intermediate subfactor 
A G C G N iV°pp where the doubled G-G system is the untwisted quantum double 
VK = K%{K) = KXk-^k{K X K). Then A C G is a braided subfactor of index 
\G/K\. The full system c^c = K%j^_^k{G x G) and chiral systems C'^c - cP^a = 
K\j^_^q{G X G). The branching coefficients and sigma-restriction K^{K) — )■ Kq{G) 
are given by {g,Ti) — ?■ T,{g' ,\Yi.dj^^^ vr), where the summation is over equivalence 
classes of g' in G conjugate to g. The trivial example K = 1 recovers the canonical 
modular invariant fll.l9p . In the case of the dihedral group G = D^, and its cyclic 
subgroup K = 7j,y, in which we are interested, the canonical modular invariant arising 
from the inclusion is given in fl3.1ip below. 



1.4 VOAs and vector- valued modular functions 

For the basic theory of VOAs, see [17]. Let V be a VOA. Write Irr(V) for the 
set of irreducible V-modules M. Among other things, V and its modules carry a 
representation of the Virasoro algebra Vir = Span{L„, G}nez, where the central term 
G acts as a scalar c = c(V) called the central charge, c also arose in Section 2.1, 
where it was defined only mod 24. The action of Lq G Vir on each V-module M 
defines a grading (by eigenvalues) on M into finite-dimensional spaces. 

Definition 3. A VOA V is called rational if 

i) V G Irr(V) and V is isomorphic to its contragredient as V-modules; 

ii) writing V = ®nezVn for the grading by Lq of V, we have V„ = for n < and Vq 
is 1- dimensional; 

iii) every weak V-module is completely reducible. 

See [17] for the details, which play no role in the following. There is no standard 
definition of rationality — we chose Definition 3 to guarantee the existence of modular 
data. All VOAs considered in this paper are rational in this sense. 
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Rational VOAs V (or if you prefer, rational CFT) realise modular data as follows. 
The primaries consist of the finitely many modules M G Irr(V) = $, and the vacuum 
is V itself. Define their characters to be the graded dimensions 

c/iM(r) = TrMg'^"-^/'\ (1.21) 

using the above grading by Lq, where as always we write q = e^'^''^. These cHm will be 
holomorphic throughout the upper half-plane H = {r G C : Im(r) > 0} [62]. Collect 
these finitely many characters into a column vector c/i(r). Then [^2] showed there is 
a representation p of SL2(Z) such that 

The matrices S, T now defined by (11. 2p constitute modular data [36]. It may or may 
not be unitary. Unlike braided subfactors, which are naturally associated to modular 
invariants, VOAs see only one chiral half of the rational CFT, and so capture only 
the notion of modular data. 

There is a rational VOA V{L) corresponding to any even positive-definite lattice L 
(the central charge equals the dimension n of L), recovering the modular data of (11.91) . 
The character corresponding to coset [v] G L* / L is c/i[t,](r) = f]{'^)~^'l2x&[v]9^^^'^- 
There is a rational VOA V{Qk) corresponding to affine algebra q^^^ at positive inte- 
gral level k, recovering the modular data of [12]. The character chx corresponding to 
highest- weight A coincides with the affine algebra character xx^ specialised to r G H. 
Affine algebra and lattice VOAs overlap for the simply-laced q at level 1, correspond- 
ing to the root lattices L. Finite group modular data V^G is recovered by taking the 
orbifold (i.e. sub VOA of fixed-points of G) of a subgroup G of automorphisms 
of a holomorphic VOA (i.e. a rational VOA V with Irr(V) = {V}) — it has been 
conjectured that any such will itself be rational, but a general proof of this still 
seems far away. Examples of holomorphic VOAs are the V{L) for L self-dual — e.g. 
for = 8 there is the Eg root lattice. In this case the orbifold theory is under better 
control, see e.g. [13], and indeed this is the main one we'll consider. A sexier example 
of holomorphic VOA though is the c = 24 Moonshine module VK 

Incidentally, Theorem 1 suggests that: a VOA corresponding to a quantum double 
should be a conformal subalgebra of a holomorphic VOA. 

Definition 4. Let p be a d-dimensional representation of SL2{'L), with diagonalT := 
pil\)- A weakly holomorphic vector- valued modular function c/i(r) with multiplier 
p is a holomorphic function EI — )■ satisfying 111.2^) . with q-expansion 

oo 

c7i(r) = cM" (1.23) 

n=0 

for some diagonal matrix A satisfying T = e^'^^'^, where each Fourier coefficient chn G 
C'' is independent of t. Let J^{p) denote the space of all such ch for p. 
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Equation fll.23p means c/i(r) is meromorphic at the cusp. [62] showed that the 
character vector of a rational VOA is a vector-valued modular function in this sense. 
These character vectors will help us study and identify the VOA. 

The characters of rational VOAs also satisfy other conditions. For instance, by 
definition fll.2ip the Fourier coefficients chn in (11.231) are nonnegative integers. The 
vacuum character c/iy = c/iq begins with Ig"'^/^^. In a unitary VOA, A in fll.23p can 
be chosen so that Xm > — c/24 for all M ^ V in Irr(V). 

Recall charge-conjugation = 3"^ from Definition 1. A consequence of fll.22p 
is that Cch = ch. For convenience we assume in the remainder of this subsection 
that C = I, i.e. we have a representation of PSL2(Z). This holds for most modular 
data considered in this paper (other than P'^Z2„+i), but at the end of Section 2.4 we 
explain how to reduce C 7^ / to the C = I case. 

Given any d-dimensional PSL2(Z)-representation p, [2], [3] describe how to find 
all c/i(r) G A^(p). In particular, there is a d x d matrix S(r) with the property 
that cA(r) G M{p) iff c7i(r) is of the form X(r) = E:(r)P( J(r)), where P{x) is 
any column vector whose entries are polynomials in x, and J(r) is the Hauptmodul 
J(r) = + 196884g + ■ ■ ■ . Thus it suffices to find S(r) for p. S(r) is determined 
through a differential equation it satisfies. This differential equation and the relevant 
initial condition is determined from two dx d matrices of numbers. One is a diagonal 
matrix A satisfying among other things the relation e^'^'^ = T; the more elusive one 
is called x- terms of A and the matrix S(r) has g-expansion 

00 00 
S(r) = 5^ q-E^ = q\lq-' + X + ^ ^"-«) ' 

n=— 1 n=l 

where S_i is the dx d identity matrix. So is the matrix x, and for n = 1, 2, 3, . . . the 
matrix S„ is recursively defined by the commutator 

[A, S„] + (n + 1)S„ = ^iifn-i{^ - /) + ^7,i-i(x + [A, x])) , (1-24) 

i=-i 

where fn,gn are defined by (J - 240)A/£'io = ^/^g" and A/Eio = Y.dnq'' for the 
discriminant form A = 77^^ and Eisenstein series Eiq = E^Eq. That is, the (z, j)-entry 
of S„ is the (i, j)-entry of (ll.24p . divided by An — Ajj + n + 1. 

The following method suffices to find A, x for the p considered in this paper 
(though invariably there are more elegant methods). First, decompose p into irreps. 
What is special about irreps p is that their A^(p) is cyclic (Theorem 4.1 of [3]): the 
space M{p) is a module over the ring C[J, Vi, V2, V3] of differentiable operators 

Vi = ^Do , V2 = ^DiDo , V3 = ^D2DiDo , for = g _ 'Ie, , 
A A A dg 

and if p is irreducible, then Ai{p) (and hence A and x) is generated over that ring 
by any nonzero ch G M{p)- All irreps occurring in this paper are subrepresentations 
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of the modular data coming from even lattices, and so the desired nonzero modular 
function can be built from lattice theta functions. 

Knowing A and x for ^ given S, T is equivalent to knowing A and x for uT for 
any third root of unity u, but the explicit equivalence [3] is not easy. Write S*^*^''(r) for 
the matrix H(r) for S,^^T and assume H'^*^^(r) (hence A*^") and is known. Then 
the columns of S^^^ are linear combinations over C of the columns of r]~^^ E'^'E^^^ and 
^-16 (^E^D^'^i^) — El'E^^^A^'^^ — 1)), while the columns of S^^^ are linear combinations 
over C of r]-^E^E'^''^ and r/-^ (DiDqS^ - E4S(o)(A(o) - J)(AW - |J)) . 

The contragredient {p^)~^ is handled similarly [3]: the columns of £'io-E'i4A^^(S-^)^ 
are linearly independent vectors in M.{{p^)^^). 

2 Comparing the Haagerup, Sym{3) and SO (13) 
2.1 The tube algebra of Sym{3) 

Recall the discussion in Section 1.2. Write S3 = Sym{3) = {e,u,u'^,T,Tu,Tu'^}. The 
three conjugacy classes are Kg = {e}, Ku = and Kr = {r, rn,rn^}. Write 

5*3 = {l,e,a} where e is sgn and a is 2-dimensional. Then each pair {Kg,TT), with 
n E Zg, gives rise to a simple component of the tube algebra. 

There are two half-braidings for endomorphism p = e as P{Ke,i) = P{K^,e) = e, 
and a third attached to p(K^,a) = e + e. We denote these half-braidings by e^^\e^'^'' 
and 2e, respectively. The conjugacy class provides an endomorphism P{^Ku,-k) = 
u + V? for every irrep vr of the centraliser — Z3, and hence three half-braidings 
(m -|- v?')^^\ {u + u^Y'^\ (u + u^Y^\ The conjugacy class has centraliser Z^ ^ Z2, 
so the endomorphism p = t + tu + tv? has two half-braidings {t + tu + tu'^)^^\ 
{t + tu + tu^Y^\ 



(c8id"")i (ugid""')! (u^^gid"")! (T8id"")i. (Tijgid-Pfjt (Ttl^gid"")! 




(l®id°PP)r (e®idPPP)i: (cr^id^PP)? 



Figure 1. Dual principal graphs for doubles of Sym{3) and Sym^S) 
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This enables us to match up Figure 1, where the bottom graph has been drawn in 
[38| Fig. 1] or [261 Fig- 31]: the upper graph arises from the Longo-Rehren inclusion 
of 5*3, and the lower describes the Longo-Rehren inclusion of the dual system 5*3, see 
[571 Remark (i) in Page 154] as well as Section 1.2 above for the general description of 
the induction-restriction graphs between A-B and B-B sectors of the Longo-Rehren 
inclusion A G B using the structure of the tube algebra. The middle vertices in 
Figure 1 also describe how the half-braidings from 5*3 and S-^ match up. 

The Longo-Rehren dual sectors associated to 5*3 and 5*3 are respectively 

[^^53] = [Ke, Indg] = 5^dim(vr)[i^e, vr] = [e^'^] + [e^^)] + 2[2e] , (2.1) 

Ks-J = Et^^^" I^d^.id] = + [1 + e] + [1 + a] , (2.2) 
9 

specialised from the discusson of Section 1.3. Observe that [^^psJ been 
computed in [381 ISSj • These braided subfactors yield the canonical modular invariants 
(recall ffTTT^ .f lT:^ ) Z22 = \cho + ch^ + 2cha\^ and Z55 = |c/io + c/ico + c/ioj^ with full 
systems S3 x S^^^, S3 x S3 respectively by [2S1 page 357], where the names of the 
primaries of VS3 are taken from the middle row of Figure 1 and the names Z22, ^55 
come from the list in [26] . 

As explained at the end of Section 1.3, the subsystem Z3 of A corresponds to an 
intermediate subfactor A G C G N ^ ]\f°pp where the dual canonical endomorphism 
of C C iV ® A^°PP is 7 = XloeZs ^ ® double C-C system is the quantum 

double = ^^^^(^3). Then A C C is a braided subfactor of index 2, with canonical 
endomorphism [6] = [0] + [b] and associated modular invariant 

Z = \cho + c/ib|^ + 2\cK\^ + 2|c/icor + 2|c/ici|^ + '^IcK^l^ . (2.3) 
2.2 The Haagerup tube algebra 

Let A = {id, u, u"^, p, pu, pv?} be the noncommutative A^-A^ system (the even vertices 
of the top graph of Figure 2) of the Haagerup subfactor A^ C M of index 5-1-1, where 
(5 = (3 -|- a/T3)/2. The fusions (product of sectors) are given by: 

[uf = [id] , [u] [p] = [p] [u]' , [p]' = [id] + [p] + [pu] + [pu'] . 

These have statistical dimensions dim(M) = 1 and dim(p) = 6. Let A = {id, a, b, c} 
be the commutative M-M system (the even vertices of the bottom graph of Figure 
2), whose fusion rules are: 

[a]' = [id] + [a] + [b] + [c], [6]2 = [id] + [c] , [c]^ = [id] + 2[a] + [6] + 2[c] , 
[a] [b] = [a] + [c] , [a] [c] = [a] + [b] + 2[c] , [b] [c] = [a] + [b] + [c] . 

Hence we obtain the statistical dimensions dim(a) = 6, dim(6) = 6 — 1, dim(c) = 5-1-1. 
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A * p pu pv? u 




A * a b c 

Figure 2. Principal graphs of the Haagerup (5 + a/T3)/2 subfactor 

The M-N sectors, [ku^] and [k,'] where k is the inclusion C M, are the odd 
vertices of both graphs in Figure 2. From this figure we see we can choose the 
endomorphism c = nuK = nu^'K. These principal graphs encode multiplication by 
K — e.g. [K][pn*] = [hi'] + [ku*]. From this we obtain the statistical dimensions 
dim(/tu*) = A and dim(K') = A (5 — 1). The remaining products M-N x A^-A^ — )■ M-N 
and M-M x M-N M-N, as well as M-N x N-M M-M and N-M x M-N N- 
N, were computed by Bisch p] and are generalised in Section 4.2 below. 

The fusions of A and A should remind one of 5*3. In particular, the system A 
has Z,3 = {id, u, u"^} as a subsystem and is some sort of perturbation of the usual 
S3 multiplication table [u]^ = [id], [u][p] = [p]^ = [id]. Similarly, the A 

system reduces to the character ring S3 when ignoring [c] , where b should be regarded 
as 5'3-representation e and a as a. These similarities were our first indication of a 
relation between the Haagerup subfactor and S3. We find in Section 4 that this 
A o 5*3 relation generalises naturally, whereas A ^ S3 does not. 

The structure of the tube algebra of the system A has been studied in [381 Sect. 
8]. Izumi introduced the following endomorphisms in S(A): 

fi = p + pu + pu^ , TTi = id + /i , 112 = 2(id) + a = u + v? + p , 

and proved that vti and 1^2 have only one half-braiding, a has three 
and finally p has six {p^^\ . . . , /x^^^) and that these half-braidings exhaust all of them, 
so that the quantum double of A has 12 primaries as in the middle row of Figure 
3. The induction-restriction graph of the Longo-Rehren inclusion of the system A is 
given in [381 Fig- 5], see also the top of Figure 3. From this we read off the Longo- 
Rehren dual sector associated to A, namely [6'x)a] = [0] + [b] -|-2[a], using the labelling 
of Figure 3. This corresponds to the canonical modular invariant (recall Section 1.3) 
^22 = Ixo + Xb + 2XaP) with corresponding full system A x A°pp, where the name Z22 
comes from the list of [27] (see also Section 3.4 below). From Figure 3 we obtain the 
quantum dimensions 1 + 35, 2 + 35, 2 + 35, 35 for b, a, c^, 0;, respectively. 
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(id8lid°P)t (iil8id°P)i (u2®id°i')i (p®id°P)i (pu(glid°P)i. (pu'lglid°P)t 




b Co Ci Cj 52 54 55 56 




(l(glid''l>)t (6»id»i')t (a^idWjt (cigid"!');: 



Figure 3. Dual principal graphs for doubles of A and A for Haagerup 

More subtle is the dual principal graph corresponding to A, on the bottom of 
Figure 3, with dual canonical endomorphism [6^^^] = [0] + [b] + [a] + [cq] and modular 
invariant Zn = |xo + Xb + Xa + XcoP corresponding to A x A°pp. To our knowledge this 
appears here for the first time. To derive it, first note that the only monomial modular 
invariants for VRg are Z22 = \cho + chb + 2cha\'^, Zu = \cho + chb + ch^ + c/ico and 
-^33 = |c/io + c/it, + 2c/zcoP (Proposition 4 below with v = 3). The number of sectors in 
A is 4, and the canonical endomorphism 7^ = (8) 0°^^ + a ® 0°^^ + 6 ® 6°^^ + c ® c°pp. 
Since (7^,7^) = 4, we have {0,0) = 4 by Frobenius reciprocity. Thus the canonical 
modular invariant Z associated to satisfies Yli^oi ~ forcing it to be Zu and 
fixing 0. Hence A has 1^ + 1^ + 1^ + 1^ = 4 irreducible sectors, call them cq, ei, 62, 63, 
one of which (say cq) is the identity Sq. To obtain the dual principal graph, in which 
Cj is connected to the primary x with (ae^yx) edges, use {ax,ay) = {0x,y), which 
holds for any type I (see Theorem 1). The fusions for the double Haagerup are 
explicitly given in [27] (or Section 3.2 below). The edges from a are determined from 
the calculations (Sa, Sq) = 6 = 1^ + 1^ + 2^ and So) = 1: call ei the sector not 
adjacent to a, and 62 the one connected to with 2 edges. The additivity of the 
statistical dimension dimSn = dim a then identifies ei = 6, 62 = 0,63 = c. Likewise, 
the edges from b come from (0:^,0:^) = 4 and (0^,00) = 1, and those from di come 
from {a^^,a^^) = 3 and (aaj,ao) = 0. From {ac-,a^.) = 5 + 6jfi, {a^.,ao) = Sj^, and 
statistical dimensions, we obtain the final edges. 

We have A C N A^°pp, the double of A, with subsystem Z3 C A. Hence by 
the Galois theory of Izumi [37], we have as in Section 1.3 and last subsection an 
intermediate subfactor A G C G N ^ ]\f°PP where the dual canonical endomorphism 
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of C C TV ® A^°PP is 7 = J2aeZ3 ^ ® ^"^^ ^^'^ double C-C system is the quantum 
double VZ3 = K'^^ilj^). Then A C C is a braided subfactor of index 2 + 35, with 
canonical endomorphism [9] = [0] + [b] and associated modular invariant (recall fl2.3p ) 

Z = \cho + c/ibP + 2\cK\^ + 2\cKf + 2\cKf + 2|c/i,j2 . (2.4) 



2.3 Haagerup modular data VHg 



The modular data VRg for the quantum double of the even part of the Haagerup 
subfactor was first computed in [3H] and simplified somewhat in [27]. It is necessary 
for all that follows though to significantly simplify it further. The result is 
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(2.5b) 



S13 



for 



for X = (13 - 3V13)/26, y = 3/V13 and c(j) = -2y cos(27rj7l3). 

It is important to note that the last 6 diagonal entries of T are To^^aj 
1 < / < 6, and the bottom-right 6x6 submatrix has entries S^i^d^, = c(//')/3 for 
1 < /, /' < 6. Our expression for that 6x6 submatrix is considerably simpler than 
the corresponding expressions in [381 1211 [35] . The proof of the equivalence is easy, 
e.g. use (11. 3p . A direct derivation is given in Theorem 5 below. 

VHg should be compared to the modular data VS^ for the (untwisted) double of 
S3, which can be obtained from (11. 8p (or see (13. ip below). As would be anticipated 
from the previous subsection, they are very similar, except that x becomes 1/2 and 
the last 6 rows/columns collapse into 2. 

The given matrix T forces central charge c to be a multiple of 24 (since Tq g — 1); 
but multiplying T by a third root of unity (and leaving 5* unchanged) allows us to 
consider the Haagerup at any multiple c of 8. Next section we generalise VUg in two 
ways: it can be twisted by a Z3 (this twist is analogous to the H^{Bq; T)-twist of 
finite group modular data, or the level k G Z>o of afiine algebra modular data); and 
it lies in an infinite sequence corresponding to the odd dihedrals. 



2.4 Clarifying the Haagerup-S'?/m(3) relation 

To emphasise that the relations between S3 and the Haagerup aren't spurious, let us 
consider possible character vectors realising these modular data (recall Section 1.4). 
This will also help identify a VOA realisation. In this subsection we focus on the 
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simplest possibility, central charge c = 8 (the smallest possible), although similar 
results hold for any other multiple of 8. 

Consider first the modular data DS^, for the quantum double of 6*3, with primaries 
0, b,a, Cj,0/ as in Figure 1. This 8-dimensional PSL2(Z)-representation decomposes 



into 3 copies of the 1-dimensional irrep with T = e 



-27ri/3 



a 2-dimensional irrep with 



kernel the principal congruence subgroup r(2), and a 3-dimensional irrep with kernel 
containing r(3). These 1-, 2-, and 3-dimensional irreps are also subrepresentations of 
the modular data for the lattices Eg, and A2®Eq respectively, and as explained at 
the end of Section 1.4 the corresponding theta functions provide all the information 
needed to extract A and x these irreps. From this we obtain for VS^ at c 
diag(2/3, 2/3, 2/3, 2/3, 0, 1/3, 2/3, 1/6) and 



=24 



the matrices A^^f 



X53 



39 


47 


81 


81 


8748 


1215 


128 


5120 


47 


39 


81 


81 


8748 


1215 


-128 


-5120 


81 


81 


167 


-81 


-8748 


-1215 








81 


81 


-81 


167 


-8748 


-1215 








3 


3 


-3 


-3 


-12 


18 








27 


27 


-27 


-27 


1458 


-152 








128 


-128 














120 


-5120 


16 


-16 














-16 


140 



(throughout the paper we write a =„ h for a = b (mod n)). 

Once we know A, x, then the simple recursion f ll.24p yields the full g-expansion of 
H(r). From H we know completely explicitly all possible weakly holomorphic vector- 
valued modular functions for this SL2(Z)-representation. Specialising to c = 8 
find there is a unique possible character vector, namely the first column of Sf|: 



we 



c/io(t) 
c/ii,(t) 
chair) = chcgir) 
ch^-^ (r) 



(1 + 39g + 699ij2 + 5761g^ + 35593?* H )\ 



v 



g2/3(47 + 671(j + 5825<?^ + 35459g 
81q2/3(i + nq + I43g2 _,_ g^y^a 
3 + 243g + 2916g2 + 21870(?=' 4 
27(?^/^(l + 22(j + 221g2 + 1476(?3 
128(j2/3(i + I6g + 136g2 + 832(?3 
I6q^/^(l + 36g + 394(?2 + 2776(?=' 



+ ...) 
+ ...) 

+ ...) 



+ ...) 



(2.6) 



) / 



It is easy to realise PS's (and hence the character vector (12. 6p ) at c = 8. Start 
with the lattice VOA V{E^) corresponding to the E^ root lattice. Its automorphism 
group is the compact Lie group £^8(IR)- Then the VOA realising (12. 6p is the orbifold 
of V{E^) by some subgroup G of -£'8(1^) isomorphic to S's (i.e. is the sub VOA of V{E^) 
fixed by G). This is indeed possible (see Theorem 4.3 of |13]). 

Now turn to the double of the Haagerup, also at c = 8, with primaries labelled as 
in Figure 3. Its 12-dimensional PSL2(Z)-representation decomposes into 1 + 1 + 3 + 7, 
where the 1- and 3-dimensional irreps are as before, and the 7-dimensional irrep occurs 
as a subrepresentation of the modular data of the 4-dimensional lattice y4352i[l, 1/4], 
in the gluing notation of Conway-Sloane [13] . From this we quickly obtain for c =24 8 
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its matrices A and x- 

diag(2/3, 2/3, 2/3, 2/3, 0, 1/3, 5/39, 20/39, 32/39, 2/39, 8/39, 11/39) 
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(2.7) 



/ 

As before, this gives us the full g-expansion of Sf^. 

There are only two possible character vectors for the Haagerup modular data at 
c = 8, namely 7 = or 7 = 1 in 



/ cho(r) \ / 

cK{t) = c/ic(,(t) 
c/ici(r) 
Chc2 (t) 
c/iai(r) 
ch„^(T) 



j2/3 
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+ ...) 



g5/39 ((7 _ ^) + (292 - 
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) / 



The proof for this (and for f l2.6p ) is similar to that sketched in Section 5.2. That 
there are sensible character vectors here is strong evidence for the existence of the 
corresponding VOA. 

We see remarkable similarities between "DHg and VS^ when c = 8. In particu- 
lar, the corresponding characters cha,chi:^ are identical, as are the sums cho + chb- 
The remaining characters c/ig^ , . . . , ch^^ of PHg can have no relation with characters 
c/io^, c/igj of T>S3, since the exponents of their g-expansions are unrelated. 

This suggests that there is a rational VOA at c = 8, call it Vs, which has conformal 
subalgebras realising the I^Hg and VS3 modular data {conformal subalgebra was 
defined at the beginning of Section 1). The Vs characters will be c/iq + cht, (the 
vacuum character) and some multiples of cha,ch^^, for the c/i's given in fl2.6p . 

Indeed, the knowledge of modular invariants says how this should go. We should 
look for modular invariants of extension type (i.e. sum of squares) for both VS^ 
and PHg, which start with |c/io + c/it,p. Both VS3 and PHg have one, namely (12. 3 p 
and (12. 4p . In the case of S3 we readily identify the corresponding chiral extension: 
Vs corresponds to the quantum double of the cyclic group Z3. The modular 
invariant says two inequivalent PZ3-primaries correspond to each cha, ch^^. 

So Vs is the Zs-orbifold (i.e. the Z3-invariant part) of the rational VOA V{Es), 
for some choice of order-3 element g of i?8(IR)- The group £'8(IR) contains four in- 
equivalent order 3 elements: orbifolding by any of these would give a VOA with 9 
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primaries realising modular data associated to a quantum double of Z3. Half the time 
this double is twisted and half the time it isn't. Take one of the order 3 orbifolds 
corresponding to the untwisted double D%j,. To recover the VSj, VOA, choose an 
order-2 element h of -^8(1^) for which gh = hg~^: the /i-orbifold of Vs realises VS^. 
Unfortunately the Haagerup VOA won't itself be an orbifold of Vs (see Section 5.1 
below), which begs Question 3 of Section 1. 

To better identify Vs, we should repeat the analysis for VZ^. Label its 9 pri- 
maries by € Z3, where (0,0) is the vacuum. Note that charge-conjugation 
C = 5*^ is no longer the identity: it sends to Hence the char- 
acters c/i(ij) and ch(^_i^^j) are equal and we should project to the subrepresenta- 
tion Span{c/i(o,o), (c/i(o,i) + c/i(o,2))/2, (c/i(i,o) + c/i(2,o))/2, (c/i(i,i) + c/i(2,2))/2, (c/i(i,2) + 
c^(2,i))/2} on which Z2 = (C) acts trivially. This decomposes into 1 + 1+3, where 
the 1- and 3-dimensional irreps are as before. Choosing the above basis, we obtain 
= diag(2/3, 2/3, 2/3, 0, 1/3) and 

(86 162 162 17496 2430 \ 

81 167 -81 -8748 -1215 \ 

81 -81 167 -8748 -1215 

3 -3 -3 -12 18 / 

27 -27 -27 1458 -152 / 

It is immediate that the only character vector possible for Vs is 

c/i(o,o) \ / g-i/='{l + 86g + 1370g2 + ll586g3 + ...) 

c/i(o,i) = ch(o,2) = c/i(i,o) = c/i(2,o) 1 ^ I g2/3(si + 13775 + 11583g2 + 71037g3 + . . . ) 

c/i(i_i) = cfe{2,2) I I 3 + 2439 + 2916(?2 + 21870?^ H 

c\i^2) = ch(2,i) / \ gi/='(27 + 594q + 5967g2 +39852q3 H 



(2i 



We recognise this as the character vector for the lattice VOA V{A2Eq), so this is the 
VOA Vs containing both the and Haagerup VOAs. This Vs also has, like V(-E's), an 
interpretation as an affine algebra VOA, and the containment Vs C V{E^) corresponds 
to the conformal embedding y42,i-E'6,i C £"8,1- It can also be realised explicitly as an 
orbifold — see |13]. The only task remaining is to identify the Haagerup VOA as a 
subalgebra of V(^2-E'6)- 



2.5 The Haagerup and SO(13) 

The relationship with concerns the primaries a, Cj and to a lesser extent 0, b. There 
also is a striking relationship with the affine algebra modular data -06,2 concerning 
the primaries 0/ (and to a lesser extent 0, b). V(-B6,2) has central charge c = 12, and 
10 primaries we'll denote 0, b = 2Ai, ai = Ag, a2 = Ai + Ag, Oi = Ai,...,95 = 
A5,t)6 = 2A6. The T-matrix is diag(— 1, —1; — — iCis )■, while the S-matrix is 
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(2.9) 
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where y = 3/vT3 and c{j) = — 2?/ cos(27rj/13) as before. Ignoring the first 4 pri- 
maries, the only difference with PHg are some signs. 

This strongly suggests relations between V(-B6,2) and the (still hypothetical) 
Haagerup VOA involving the Goddard-Kent-Olive coset construction [33]. In the 
VGA language, the coset construction was developed in |29]; see also the lucid treat- 
ment in Section 3.11 of There are several ways this could go: e.g. the Haagerup 
VOA at c = 8 could be a coset of V(-B6,2) by a c = 4 sub VOA V4. In this case the 
characters of V(-B6,2) (which can be determined by the Weyl-Kac character formula 
[I2I Ch.l3]) would be built from the characters of V4 and the Haagerup. 

Although the possible relation between the Haagerup and V(-B6,2), in particular 
those involving the coset construction, seems very intriguing, for reasons of space we 
limit the discussion in this paper to relations between the Haagerup and S3. 

3 Generalising the modular data of the Haagerup 

3.1 Dihedral groups and orthogonal algebras 

Last section we observed that the Haagerup modular data is closely related to that 
of the symmetric group 5*3 = Sym{3) and afiine so(13) at level 2. More generally, we 
propose next subsection a two-parameter generalisation of the Haagerup, related to 
the odd dihedral groups and affine so(2m + 1) at level 2. 

First, we compute the modular data for quantum doubles of the dihedral group 
Di, = (r, M I t'^,u'^,tu = u~^t) where u = 2n + 1. Of course 5*3 = D^. The twist 
group is H^{Bj:,^;T) = l^^v — ^2 x Zj,; because the Schur multipliers of Dy and 
cyclic groups all vanish, the modular data is cohomologically trivial and is given by 
(11. 8p . The conjugacy classes of have representatives e,u^,T for 1 < h < n, with 
centralisers D^, (u), (r) respectively. There are two 1-dimensional irreps of D^, (call 
them eo = l,ei) and n 2-dimensional ones (call them cTj for 1 < i < n); denote the 
1/ 1-dimensional irreps of (u) = by ttj for < j < u, and the two 1-dimensional 
irreps of (r) = Z2 by eo, ei again. The primaries fall into four classes: 

1. two primaries: the vacuum := (e, 1) and b := (e,ei); 

2. n primaries, labeled ai := (e, (Jj) for 1 < i < n; 

3. nv primaries, labeled C/ij := (u^, vTj) for 1 < /i < ri, < j < i/; 
Ji^. two primaries, labeled Hi := (t, e;), I = 1,2. 

Fix a sign s = ± and integer cu G Z. Corresponding to these 4 classes we have 
the modular data 



T=diag(l,l;l,. 



. . . , 1; exp[27ri (w/i^ + z//ij)/i/^]; t, —t) , 



/ 22x2 
^ lnx2 2. 
I' lnvx2 



\ pB^ 




(3.1) 
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where t = l,i for s = 1, — 1 resp., kaxb for any number k is the a x b matrix with 
constant entry k, 

Di,{h,j) = 2 cos{2tt ih/u), and E(^hj)^(^hij') = 2 cos(27r(2a;/i/i' + uhj' + vh'i)/v'^). We'll 
denote this modular data by V^'^D^. The untwisted double is V~^'^Di,. Note that 
'ps,uj+kujj^ is equivalent for any integer k, so the twist does indeed live in Z2u- We've 
separated the order 2 from the order u twists, for later convenience. The quantum- 
dimensions (=statistical dimensions) S'*o/5'oo of b, Oj, Chj, f « are 1, 2, 2, z/, respectively, 
and the global dimension 1 / 6*0,0 is 2z/. 

The S,T entries live in Q[.^4i/2]. Choose any Galois automorphism a G 
Gal(Q[,^4^2]/Q) fixing 1 and define i G by a^^2 = ,^^2- Then applying a to 
T>^'^Dy, as explained in Section 1.1, yields modular data equivalent to V^'^'^D^. In 
particular the fusion ring of 'D^^'^Dy depends only on z/, s and gcd(a;, v). In fact, if 
£ G Z^, then V'^^'^D^ and V'^^D^ are equivalent. 

The second ingredient going into our generalisation of the Haagerup modular 
data is the affine algebra modular data -Bm,2; which has central charge c = 2m. In 
particular [H], -6^,2 has m + 4 primaries, labelled 0, b = 2Ai, ai = A^, 02 = Ai+A^, 
di = Ai for 1 < I < m and dm = 2 Am- Write fi = 2m + 1. The T-matrix is 
eiTdiag(l, 1; -il^- if), while the ^-matrix is 

/ 22x2 

S= \ 

where x = 1/y^, B and F are as above, and Hi ir = 2xcos(27r//'//i). The quantum- 
dimensions of b, ai, Di are 1, ^/JI, 2, respectively, with global dimension 2y/fl. 

There is no obvious twist of the -Bm,2- Curiously [15], affine so(z/^) at level 2 
coincides with V^D^ for a specifically chosen twist u (this is one of the few overlaps 
of affine algebra and finite group modular data [ISj), but that fact seems to have no 
role in our story (in particular we are not interested in the values m = [u"^ — l)/2). 

3.2 Generalising the Haagerup modular data 

Sections 2.3, 2.5 and 3.1 suggest a generalisation of l^Hg. Write u = 2n + l for n > 0, 
and choose G Z as before. Write m = 2{rP + n + 1), /i = 2m + 1 = z/^ + 4, and 
S = + y/Ji)/2. Note that 6 satisfies 6"^ = uS + 1 and lies in the range u < 6 < u + ^. 
Again the primaries fall into four classes: 

1. two primaries, denoted and b; 

2. n primaries, denoted Oj for 1 < i < n; 

3. nu primaries, denoted Chj for 1 < h < n,0 < j < u; 
4- m primaries, denoted 0/ for 1 < I < m. 



F 02xm, 
0mx2 H 



(3.2) 
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Breaking S and T into 16 blocks, as in the previous subsection, we get 

/,,2i. „,,„ro„, ;2_ / 



dia 


g(i,i;i 


. . . , 1; exppTTi 






( ^ 




B' \ 
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lnx2 


2nxn 


Onxm 
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-pH ) 



(3.3) 



where fcaxft, -D, -E and are as in Section 3.1 (so —uHi^i' = — 2y cos(27r//'/yu)), and 

, 1 fl-y l + y\ , f 1 1 ■■■ 1 

A = - i ^ . and 5' = w I 

2 Vl + y 1- V V-1 -1 ■■■ -1 

for y = Denote this modular data by T)'^B.g2n+i', we call V^Hg^ Haagerup- 

Izumi modular data. To our knowledge this modular data is new except when n <1 
and V divides in particular twisting the Haagerup data was unanticipated in the 
literature. We discuss further generalisations next subsection. 

The quantum dimensions (or statistical dimensions) S'*,o/'S'o,o for Oj, C/ij, 9/ are 
respectively l + z^(5, 2 + i/(5, 2 + z/5, and p5. The global dimension is 1/5*0,0 = ^ {yd + 2). 
Note that the submatrix —vH has exactly half of the // — 1 rows and columns of S. 

The fusions can be computed directly using Verlinde's formula ( 11. ip (or from 
Proposition 2 below). As mentioned in Section 1.1, it suffices to consider all /c 
different from 0. For P'^Hg^,, all of those fusion coefficients are 1 except for 

iVb,ai,a, = 2 , iVb,Ch,i,c;,,, = 2 , A^b,0,,0, = , iVa„a,,afc = 2 H k =^ SI + sj , 

= 2 if h" =u sh + s'h' , k=usi + s'j + 2u {sh + s'h' - h)/p , 
^a„ch,,,c;,,fe = 2 if sz =^ j - A; , iV5,,o,,,o,„ = if /" =^ s/ + s'l' , 

where s, s' G {±1} are arbitrary signs. 

As in Section 3.1, 'D'^}lg2n+i depends up to equivalence on the value of the twist 
uj mod z/; the fusion ring depends up to equivalence only on n and gcd(i^, w); and 
V^Y{g2n+i — '^Hg2n+i ^uy i coprime to z/. 

The computation of fusions reduces to the identity 

A / 2Trad \ ( 2iTbd \ ( 2Txcd \ , 
8 > cos — cos — cos — = -4 + (2A; + l)s 



proved using 2 cos(x) = e'^ + e where < s < 4 is the number of pairs (s^, of 
signs ±1 such that c =2fc+i s^a + s^fe. Modularity reduces to the Gauss sum 

2fc 

for k even 



^exp[27ricaV(2A; + l)] = v/2^(^^^^ || 



for k odd 



Since /i =8 5, /i can never be a perfect square, and the Jacobi symbol j will always 
equal —1. See Proposition 2 below for a more elegant argument. 
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Of course P^Hgg recovers the original Haagerup double T>}ig, given in Section 2.3. 
The similarity with the modular data V^'^^Dy and 5^,2 is evident (ignoring the first 
4 primaries of -Bm,2; the only difference with V^B.g^, is that the bottom-right corner 
of both T and S are off by some 6th root of 1). In particular, class 1 runs through the 
2-dimensional irreps of Di,, while h and j in class 3 parametrise the size-2 conjugacy 
classes in D^, and the irreps of Z,y, respectively; class 4 runs through the fundamental 
nonspinorial weights of so(/i). 

The mysterious '13' in "DHg is thus 4 + 3^, where the 3 here references the normal 
subgroup Z3 of 5*3. It is tempting to guess that the 4 is 2^, where the '2' is the 
involution of S3. This suggests the further generalisation of P'^Hg^ given in the 
Wildly Optimistic Guess of Section 1. 



3.3 Further generalisations 

The way in which P'^Hg^, is built from V~^'^Dy and Bm,2 leads to the notion of grafting 
modular data. A simple instance is provided by the following proposition, but it can 
be massively generalised. For instance it would be interesting to extend it to vacuum 
blocks larger than 2x2. 

For simplicity restrict here to unitary modular data (recall Section 1.1) with real 
matrix S. Call modular data ($,0, S*, T) Z2-/ammate(i if it can be written in block 



b,b 



To,o and 



form $ = (0, b; ai, . . . , a™; ^i, . . . , Dn) such that S'o,b = 5*0,0: 
Sb,(3j < 0. It is elementary to verify that modular data is Z2-laminated iff it can be 
written in the form 



^ X 


X 


a 


d \ 


X 


X 


a 


-d 


-'T 

a 


(f 


A 


Omxn 






Ojnxn 


D J 



T = diag(r, r; s; t) 



(3.4) 



for some numbers x, r G M, row vectors a, G 



d,t e 



and matrices A, D. 



The condition S^fl = Sq^ implies the matrix Ni, with entries A^J ■ for all i,j & ^ is 



a permutation matrix (see e.g. [31]); in conformal field theory, such a primary b is 
called a simple- current. The modular data V^'^D^ and -8^,2 are both Z2-laminated, 
for any odd v and any twist (±, uj) and any rank r, as is Ai ^. 

Proposition 2. Consider Z2-laminated modular data ($, 0, S, T) and ($', 0', S', T'). 



(a) There are integers M>0, L>0, A^>0 such that 2x = 1/y/M, N and 4^ divide 
AM, 2^^^^/Md G Z"", 2^/W/Nde Z", and the gcd over all components 2^~^y/Mai 
is 1 (likewise for the components 2a/ M/Ndi^). 

(b) Suppose X = x' (that is, the global dimensions agree) and r = r' (i.e. the central 
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charges c, c' agree mod 24)- Then the following defines modular data: 



S 



S 





^ X 


X 


a 








X 


X 


a 


-d' 










A 









yd"" 







D') 


and r 


= —r'u for some 


Then 










fx. 


x+ 


a 


a' \ 






X- 


a 


-a' 








A 














-A') 





T = diag(r, r; s; t') . 



T = diag(r, r; s; —ut') 



satisfies all conditions of modular data except possibly that the fusion coefficients are 
nonnegative integers. Its fusion coefficients are all integers iff 4\{M' — M), and both 
L > 0, L' > 0. S,f define modular data if in addition N', , , < 8^'/(M' - M). 

Proof. Recall (11 ■4p . Clearly, any a permutes {0, b}, so the number will be 
fixed by any a. But is an algebraic integer, being the sum of squares of certain 
eigenvalues 5'j^o/'S'o.o of integer matrices. Therefore is indeed an integer. Com- 
puting A''o„o^,B;^/A^Q^.,E,;,,afc = ai/aj and iVa„Ofc,Ofe,/^'^n,,£)fc,Ofc = we see that there are 
/3, 7 > for which both vectors /3a and 7(i are integral. The gcd condition follows 
by choosing these /3,7 as small as possible. Because d ■ d = 1/2, we know 7^ G Q. 
Because Na^^Of^^i,^ = 2aid\/x G Z and each quantum-dimension dk/x is an algebraic 
integer, we know ai/x G Z and hence {fix)~^ G Z. From 2x^ + a- a = 1/2 we now see 
that 4 divides and (/3a;) ~^ is a power 2^ of 2. This gives us (a). 

The proof of (b) and (c) is now straightforward. Most interesting are the fusion 
calculations in (c). Note that Ax^x''^/{x^ - x'^) = 1/{M' - M), ai/x G 2^Z and 
a'l^/x' G 2^'Z. We obtain 



b,b,b 



M' 



2a, 



N, 



2al 



X' {M' - M) 



M 



b,b,ai 



X {M' - M) 



xx' {M' - M) 

aiiCLjiCii" 



N, 



x2 (M' - M) 
x'2 (M' - M) 



Si 



-6, 



k,k' 



2x3 (M' - M) 
2xx'^ {M' - M) ' 



2x2^' (M' - M) 



(Xp,(ly(lp,n 



2x'3 (M' - M) 



QED 

We call S", T the gra/i of S' ,T' onto S", T. Of course V^Yig^ is the graft of 5^,2 
onto V^'^D^; here M' = /i, M = i/^ (so M' - M = 4), w = +"-^ and L = L' = 2. 
An easy example of (b) is grafting V^'^Di, onto T?^''^Di, to form 'D~''^Dy. 
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An easy generalisation of V^Hg^, is to replace with any abelian group K of odd 
order u. Take a trivial twist a; = 0, for simplicity. The entries of S and T involving 
classes 1 and 4 are identical to those in f l3.3p . The role of D^, is played by the semi- 
direct product Dx = /'^xiZ2, so the class 2 primaries are labelled by the 2-dimensional 
irreps i of Dk, and the class 3 ones by pairs {h, j) where h runs through representatives 
of the cardinality-2 conjugacy classes of Dk and j runs through the irreps of K. The 
T-entries for class 2 are again I's, while those for class 3 are the evaluation j{h). 
The submatrices D and E are again cosines, namely Di^(^hj) = i{h) and E(^hj-^^(^h>j/-j = 
2Ke{j'{h) j{h')). Although the modular data thus generalises naturally to arbitrary 
odd order abelian groups, we see in Section 4.1 that the subfactor realisation for 
noncyclic K is more subtle. 

It is tempting to look for additional twists of VRg^^. After all, we've suggested 
that the generalised Haagerup is a twin of the dihedral -Dj,, and the latter can be 
twisted by Z,2iy and not merely Zj,. The Z^, twist of both VDi, and T>B.gj^ affects the 
class 3 primaries. The independent Z2 twist of VD^, affects the class 4 primaries, so 
this suggests one should look for additional twists (by Z2 or perhaps Z^) of VB.gj^, 
which sees the m primaries of class 4 and leaves the other classes untouched. 

We haven't gotten this to work. For instance, one approach to probing some of 
the missing twists would be a Galois automorphism acting on S and T entry-by-entry. 
The entries of V^B.g^ lie in Q[^i/2^], with Galois group Z^2 x Z^. The class 3 twists 
uj & are stable under Z^a, as explained last subsection, so we should consider the 
effect of on P'^Hg^. 

For i E li^, aey = ±y for some sign. If agy = —y, then ag sends 'D'^Hg^ to 
nonunitary modular data with exactly the same fusion coefficients: the vacuum pri- 
mary is still the first one, but the positive column of S is the second one. This means 
this modular data won't have a subfactor realisation (it could perhaps have a planar 
algebra interpretation, if the requirement of positive-definiteness there is dropped), 
though it can still have a VOA one. The notion of twist should presumably pre- 
serve unitarity, so we won't regard this aiV^Hg^ as a twist. Nevertheless it may be 
interesting to search for its VOA realisations. 

For this reason we should focus on those £ G Z^ for which aey = +y. Now, it is 
typical for small n that is a prime power (this is true for all u < 19 except u = 9). 
In this case, any such ae will send V^B.g^^ to equivalent modular data. The search 
for new twists using Galois seems to fail. 

However when fi is not a prime power, this idea bears fruit. Given the prime 
decomposition /i = nLi^rS where each > 0, then a£ E Gal(Q[,^,^2^]/Q) = Z^a^ 

fixes 5 iff Y[i=i (^^j — +1- We can also require here that i =^2 1 (the z/^-part of 
£ merely shuffles the twists u). The automorphism an maps 'D'^Y{g^, to equivalent 
modular data iff = +1 Vz. Given any set P (possibly empty) containing an 

even number of distinct prime divisors of /i, pick any Ip G ^^2^ with £p =,^2 1 and 

= —1 iff Pi G P. By P'^'-^Hgj^ we mean the modular data aipS,aipT; up to 
equivalence it is well-defined, and is inequivalent to P^^Hg^ even though it has the 
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same fusions. Thus as long as > 1, i.e. /i is not a prime power, this construction 
yields new modular data. In Section 4.1 we find subfactor realisations for this new 
Galois twist, at least when u = 9 and u = 19. 

We learned in Section 2.4 that l^Hg also sees VZ^, in some ways more directly 
than it does VS^. Since VZ^, has only a Z^^ worth of twists, it is certainly not 
inconceivable that the correct generic twist group for "PHg^ is indeed Z;,. 

3.4 Miscellanea involving modular invariants 

[26| Sect. 6] studied the quantum double S3 modular data in detail, and similar 
techniques [28] can be used for their twists, as noted in [21]. Using [271 Theorem 
4.3], [28] (with the correction to [27] where Zn was erroneously recorded as being 
sufferable, when we have now seen in Section 2.2 of this paper that it is indeed a 
dual canonical modular invariant), we can summarize these results in the following 
table [nimless means there is no compatible nimrep, while nimble means there is; 
sufferable means there is a subfactor realisation): 





nimless 


nimble, insufferable 


sufferable 


total 




14 


6 


28 


48 




1 





4 


5 




1 





8 


9 




9 


7 


12 


28 




? 


? 


8 <? < 13 
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Table 1. Modular invariants for the doubles of Sym{3) and Haagerup 

The commutant of S, T of course forms an algebra under matrix multiplication, 
but the sufferable modular invariants themselves have a fusion structure, in the sense 
that the product of two of them is always a linear combination over nonnegative 
integers of other sufferable modular invariants. The modular invariants for V^'^D^ 
and V^'^'^Dy, including this fusion structure, are naturally isomorphic for any £ G Z^, 
using aiV'^^D^ ^ V'^'^D^; the same applies to V^Rg^ and V^'^Rg^. 

[28| [2Tj remarked that there is an (injective) homomorphism from the 28 "PHg 
modular invariants to the 28 sufferable V'^'^Ss ones. We can define a natural injective 
homomorphism 4>^^gl'^^^ from the sufferable twist s,uj quantum double S3 modular 
invariants into the untwisted 5*3 ones. We can actually embed the sufferable —1 into 
+1 modular invariants, the sufferable —1 into the —0 modular invariants and the 
sufferable —0 into +0 modular invariants. 

This can be generalised and understood as follows. As is clear from Sections 
3.1 and 3.2, primaries for both V^^'^D^, and D'^Hg^ fall naturally into 4 classes, and 
so their modular invariants decompose naturally into 4x4 blocks as did their S, T 
matrices. Write 5'[j]^[j], Tjjjj,], M[j]jj] for these blocks. 
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Theorem 3. Write k for the number of distinct prime divisors of jj,. Define a map 



( M[i]ji] M[i]j2] Af[i]j3] 

^[2],[1] ^[2], [2] ^[2], [3] 

^[3],[1] ^[3], [2] ^[3], [3] 

V 02x2 02x„, O2 



(a-6)/2x2 / 




where 



a b 
c d 



) 



Then ^(p- 



is a bisection between all modular invariants 



of X^'^Hgy with b ^ 0, and all modular invariants of V'^'^Dy with b ^ 0. "^^Hg!^ 
is a 2^~^-to-l surjection from all modular invariants of V^B.g^^ with b = 0, onto 
all modular invariants ofV^^'^Dy with 6 = 0. Moreover, any modular invariant of 
V^'^^Dy is a modular invariant ofD'^^'^Dy (but not conversely). 

Proof. Let M and M' be modular invariants for V^Hg^^ and 'D^^'^Dy, respectively. 
Let 5, T (resp. S',T') denote the modular data for V^Hg^ (resp. V^'^Dy). From 
M'T = T'M' we see that Mj'.j = = Mf^j [.] for i = 1, 2, 3. We also see that M^'^j 
is diagonal. Recall fll.l2p : we see from the 0-row of S that in classes [1],[2],[3] is 
identically +1 for all £ G ^^2^, while s^{^i) = j equals —1 when e.g. i =^ m. Thus 
also M[j]j4] = = M[4]^[j] for z = 1, 2, 3. 

Evaluating MS = SM at ([1], [4]) and ([4], [1]), we get two possibilities (identical 
conclusions hold for M'): 

(i) M[i].[i] = / and M[4]j4] is a permutation matrix: Mg^^o^, = ^jr/,/' for some permuta- 
tion vr of 1 < / < m; 

(ii) M[i]ji] = l2x2 and M[4]j4] = 0. 

Looking at the remaining equations MS = SM, MT = TM, M'S' = S'M' and 
M'T' = T'M', we see that (in case (i)) the equations involving M[4] [4] decouple from 
the others and reduce to 



On the other hand, the equations for the unknown entries of M[j] yj, i,j 7^ 4, are iden- 
tical to the corresponding equations for Mj'^j yy For example, (MS') [2]j4] = (S'M) [2]j4] 
and (M'S")[2],[4] = (S'M')[2],[4] both reduce to M„^,o = and M^^ g = M'^^ ^,. There- 
fore in case (ii), where M[4] [4] = = Mj'^j j^j, this means "^^Hgi^ is indeed a bijection. 

In case (i), we know M'^j = I, so we need to identify the permutation vr. Write 
ttI = i. Then 03.51) requires i"^ =^ 1, hence nl =^ So to any 1 < i < m satisfying 
£^ =^ 1, define nil to be the unique number 1 < nil < m obeying tt^I =^ ±li. Then ni 
satisfies (13.51) and thus defines the M[4] [4]-block of a case (i) V^Hg^ modular invariant. 
These 2^~^ i parametrise the kernel of '^4>]lgy in case (i). 

V+'^^Dy is handled identically. QED 

Note that ^(p^gy is linear and preserves matrix multiplication. In particular, when 
/i is a prime-power (as it is for the original Haagerup PHg, and all < 19 except u = 



7i{l)n{l') =,±ll', Tx{lf =,P. 



(3.5) 
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9), there will be an algebra isomorphism between the span of modular invariants for 
P'^Hgj^ and for V~~'^Di,, and an (injective but nonsurjective) algebra homomorphism 
of these into V+'^^D^. 

The relation between the 28 sufferable V'^^'^S-s modular invariants and the 28 
P^Hgg modular invariants would seem to be a coincidence: e.g. at z/ = 1, V~^'^Di = 
P+''^Z2 has exactly 6 modular invariants, all sufferable, while V^Hg^ has exactly 2. 

The dihedral group is only half the story; the other half is the affine algebra data 
i?m,2- Its modular invariants were classified in [32]: when /x is not a perfect square (the 
situation here), the complete list is B{do,io) and B{di, ii, d2, ^2) = B{d2, ^2', di, ii) 
where di\fi, fi\df, and 1 < £i < m obeys i'f =ifi/^ 1- for matrices B{d,£) and 
B{di,ii]d2,i2) defined in [32]. Those modular invariants which are permutation 
matrices are the B{^,i) for 1 < i < m satisfying i"^ =^ 1: their nonzero entries are 
B{fi,i)z,z = 1 for 2; G {0, b, tti}, and B{^, £)^^^j,^, = 1 for /' =^ ±£l. We see from this 
and Theorem 3 that the modular invariants of V^B.g^ come from those of V^'^Hg^^ 
and Bm,2 in a very direct sense. 

Recall [531 122] that the sufferable modular invariants of the double of a finite 
group G are parametrised by pairs [H, ip) where if is a subgroup of G x G and 
ip G H'^{Bh]C^) is called discrete torsion. This suggests: 

Question 4. Find some analogue for P°Hg^ of this {H, ip) parametrisation of suf- 
ferable modular invariants. 

This {H, ip) parametrisation belongs most naturally to the VA formulation of the 
double of the group G. As mentioned in Section 1.2, VA is also the double of G. Is 
there also a parametrisation of the sufferable modular invariants of the double of G 
which is more natural in the VA language? 

Of special importance (recall Section 1.3) are the monomial modular invariants 
of V^B.gj^. Three of these are obvious: 

|c/io + c/if, + 2^c/i„,p, |c/io + c/i(, + ^c/ia^ + ^c/ich |c/to + c/tb + 2y^c/ic;,op. 

(3.6) 

Proposition 4. When u = p orpq for (not necessarily distinct) primes p, q, the only 
monomial modular invariants are the three in (13.61) . 

Proof. We need to find all eigenvectors u of S,T with eigenvalue 1, with uq = 1 and 
all other G {0, 1,2,.. .}, so Z = | u^chxl'^. We know = 1 since M(,^(, = 1 for 
any modular invariant (clear from the proof of Theorem 3), hence Mg, = for all /. 
We need to determine Ui := and Uhj := u^^ 

Tu = u implies Uh,j = unless v divides hj. Therefore if Uh,j 7^ we have: 

(i) for V prime, j must be 0; 

(ii) for V = p"^, only for j = or p\h and p\j] 

(iii) for v = pq (p 7^ g), only for j = 0, or p\h and or q\h and p\j. 
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Su = u implies, for all 1 < h, h" < n and < i < 

2 + 2 Uj + 2 Uh'j cos(27ri/i'/ z/) = z/Wj , (3.7) 
2 + 2^cos(27rj/i/z/) + 2 ^ cos(27r(/ij + h'i)/u)uh'j = h'Uh,i , (3.8) 
cos(27r/ij/ z/) M/j// J = cos{27ih"j/ v) Uhj , (3.9) 

where we sum over h',j (f l3.9p was obtained by hitting (13.81) with cos(27r/i"i/z/) and 
summing over i). In analysing these equations, it is useful to recall (for k odd) 

= 2 5^cos(27r/A) = /i(A:), (3.10) 

the sums taken over the 1 < i < k and 1 < / < A;/2 which are coprime to k. The 
Mobius function fi{k) equals (—1)" if A; is a product of n > distinct primes, and 
otherwise. The Galois symmetry ae for £ G (recall (11.111) ) maps to a±ei and Chj 
to c±ih,ej. Therefore Ui = Ugcd{i,u) and Uh,j = Ugcd{h,u),jh/gcd{h,u), so Uh-j = Uh,j. 

For u prime, this means Ui = Ui and Uh,j = Mi,o^i,o- Plugging this into (13. 7p gives 
Ml + Mi,o = 2, which correspond to the three solutions (13. 6p . 

For V = p^, we read off from (13. 9p with h = p, h' = 1 that Mpjp = Mp,o — Wi^o for 
all j 7^ 0. Then comparing (13.71) for i = l,p gives 2 = Up + Ui^ = Ui + Up^. We can 
now force Up = Ui by (13.81) at h = p,i = 1, and again we recover (13.61) . 

Finally, turn to z^ = pq. (13.90 with h = 1 and h' = p,q gives respectively Up = 
^ip,o - ^ii,o, Ug^pj = Ugfi - uifl for j ^ 0. Then ([321) at z = l,p, g gives + Uq,o = 
Up + Up = 2 and Ui o = 2 + Ui — Up — u^. Now (13. 8p at h = i = p and h = i = q give 
Ml = Mp = and Mi,o = Wp,o = Wg,o and we're done. QED 

The z/ = 3 modular invariants (I2.3p . (l2.4p generalise, for any V^'^Dy, V^Hg^^, to 

\cho + c/ibP + 2 ^ \cKf + 2 5^ r . (3.11) 

This corresponds to the subsystem in both VD,y and X'^Hgj^, and to the VOA 
realising V^Z^ and containing both V'^'^Dy and I^'^Hg^ VOAs. 

4 Subfactors for Haagerup— Izumi modular data 

4.1 Izumi's subfactors and their modular data 

In Section 5 we address the question of realising this generalised Haagerup modular 
data 'D'^Hgj, by VOAs. In this section we address their subfactor realisation. 

In Section 7 of [SB], Izumi suggests associating subfactors to any odd abelian 
group K of order z/, using endomorphisms in the Cuntz algebra O^^i. (Warning: 
Izumi's n is our v = 2n + 1.) Let A{K) be the set of all u x u complex matrices 
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A = (Ag^h), g,he K, satisfying 



Ag,h = Ah,g , Ag^h = A_h,g-h = Ah-g-g , Agfl = Sgfi - 1/(5 - 1) , (4.1) 

Ag+m,hAm,h = Sg^o - Sh,o/S , (4.2) 



Am,g+hAm+k,gAm+l,h — Ag+l,kAh+k,l — ^g,0^h,o/ ^ , (4.3) 



for all g, h, k,l E K (recall 6 = {u + y/u'^ + 4)/2). These equations imply 

\Ag^h\ = V6, (4.4) 

Ag + k,hAh,k = A_k,gAg^h-k , (4.5) 

for any g,h,k G K with g ^ 0, h ^ 0, g ^ h. Izumi shows that to any matrix 
A G A{K), there corresponds a (nonbraided) subfactor of index 6 + 1 with principal 
graph the z/-star (see the top graphs of Figures 2 and 4). We call these subf actors of 
Izumi type K . Next subsection we determine much of the data of these subfactors. 
His Theorem 8.4 identifies the modular data S", T of the even part of the quantum 
double of this subfactor in terms of m(z/^ + 1) variables ui, C\,h^ g,h E K, 1 < I < m, 
satisfying the m(z/^ + 1) equations: 

^0,9 = ^' - ' ^'^i^^ - '^9+k,2hClk = Sh,o^i/S ■ (4.6) 

geK k£K 

Some solutions to 04. 6p . occurring when = 1, are redundant (i.e. correspond to 
b, tti or C/jj) and should be ignored. In particular, precisely o"i(z/) — 1 (i.e. the sum 
of divisors d > 1 of i^) solutions to (14. 6 p with ui = 1 are redundant, and for any 
other root of unity ui of order d dividing u, the number of redundant solutions is 
'Y^d' (f){v / d') /2 where is Euler's totient and the sum is over all divisors d' < d oi d. 
The entries of the S and T matrices equal those of V^Rg^ except possibly for those 
in the bottom-right m x m block, which are given by Ta,,o, = e-2--/24a;; and 



L- Lu, + S Y: CL^^C'-.^ ■ (4-7) 
\ g,heK / 



When some ui has more than one solution in (14. 6p . (14. 7p is ambiguous as it isn't 
obvious which solutions to (14. 6 p to use. 

Izumi shows in his Example 7.1 that ACZ^) contains exactly two matrices, both 
corresponding to the Haagerup subfactor. In Appendix C he exactly solves (14. 6 p for 
A G ^(Za), and in this way obtains (a complicated expression for) VHg. 

In his Example 7.2, Izumi shows A{'L^) contains exactly 4 matrices, again cor- 
responding to a single (new) subfactor. Izumi does not solve (14. 6 P for them. This 
analysis of (I4.ip - (l4.3p for z/ = 3 and z/ = 5 is as far as Izumi went (this is what we 
mean when we call it a hypothetical family of subfactors). 
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Nevertheless we can push a httle further this analysis. The group Aut(i^') acts 
naturally on the matrix A by shuffling the entries: a G Ant{K) sends A to A'^ defined 
by Ag j^ = Aag^ah- For example, when K = Z^, Aut{K) = acts by multiplication. 
In this way, given any A G A{K), Aut(i^') embeds into Aut(Q[/l]) and fixes S, where 
Q[A] is the field generated over Q by all entries Ag^h- Call A, A' G A{K) equivalent if 
A' = A°' for some a G Aut(-R'). Equivalent matrices give rise to equivalent subfactors 
and equivalent modular data. For later convenience, for x G Q[A], write Trx for the 
orbit sum ^j^g^^,) y over the -fT-orbit (x). 

It is very difflcult to solve (14. 6 p directly (see Appendix C of [38] for the u = 3 
argument, which is already quite involved). Our strategy for finding the modular 
data is to reduce it to computer calculations. A key observation is that if S, T are 
the modular data of P°Hg^, and if S',T is a second modular data, with 



q/ _ q 



1 . ^ 

< 4.8 

18^ ^ ' 



for all 1 </,/'< m, and all other entries of S' and S are equal, then (II. 3p implies 
S' = S everywhere. Maple 11 (x86 64 Linux) was used in these calculations. It 
is difficult to rigourously analyse the resulting error, but various consistency checks 
(e.g. unitarity of the numerically obtained S' matrix, or the values of l^j^ g^, — 50,^0^,1) 
all indicate it is on the order of 10~^ or better, far more precision than is actually 
needed. More details are given below. Computer calculations were not used in the 
determination of A{K) in Theorem 5. 

Theorem 5. (a) There is a unique A G -4.(1). Its modular data is I^^Hg;^. 

(b) There are precisely 2 matrices in ^(Zs). They are equivalent and realise V^Hg^. 

(c) There are precisely 4 matrices in A{1,^). They are equivalent and realise P^Hgg. 

(d) There are precisely 6 matrices in A{Tjj). They are equivalent and realise "D^Hgy. 

(e) There are precisely 12 matrices in A{'Lq). They form 2 equivalence classes and 
define inequivalent subfactors; one realises the T matrix of P'^Hgg, and the other 
realises the T of 'D'^'^^'^'^^l^gQ (defined in Section 3.3). 

(f) A{Z3 X Z3) is empty. 

Proof. Let's begin by identifying A{K) for these groups (done in [38] for K = Z3, Z5). 
It is trivial that ^(1) = {(fEf)}; ^iid for Z3 and Z5 (14. ip fixes everything in terms 
of complex numbers a and a, b, respectively: 



5- 2 -1 -1 \ 1 
^(Zs) = -1 -^1 J , A{Z,) = — 



/ 5-2 -1 -1 -1 -1 \ 

— 1 —1 aba 

-1 a -1 b b 

. 1 b b -1 a . 

y -1 a b a -1 J 



For z/ = 3, (14. 3 p forces a"^ — a + 6 = 0. The two roots of this polynomial are complex 
conjugates, and indeed — 1 G Z3 here acts as a O a. Similarly, the generator 2 G Z5 
acts on A^Z^) by sending aH-^H-ai— i-fei— j-a. Then a has minimal polynomial 
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— 6x^ + (35 — 2)x^ — 6^x + 6^ = and the other 3 roots of this polynomial are 
a,b,b — indeed, the Galois group over Q[6] of this polynomial is with this action. 
Given a, this and the equation 1 — a — a + ab + ab = then determine b. 

Computation of Ai^Zj) is handled similarly. (14. ip and (14. 5 p give us 



A{Z7 



S-1 



( 



5-2 
-1 
-1 
-1 
-1 
-1 
-1 



-1 
-1 
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b 
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b 

a 



-1 
b 
d 
d 
b 



A/ 



(4.9) 



The generator 3 G sends a\-^c^d\-^a^c\-^d^a and b ^ b ^ b. (14.5 
gives ad = be, while (g^D with {g, h) = (1, 1), (3, 1), (2, 3) give 

(6 + 6 + 2)(c + c) = -1 + Tra = {b + b + 2){a + a) = {b + b + 2){d + d) , 



where the last two equalities come from the Z^ symmetry. 

,3 



a" and we get two incompatible equations for s : 



(5 



^ -2, then 
a + a, namely 
3)£-5s2 = 55 + 1. 

and Tr a = 1. Let t = a + d + c. Then 
l + (6 + l)(3-5) and we 



„ If ^ + ^ 

a = d = c and Sb = 

+ s^{5 - 4) - 5s = 5 - 2 and 
Thus b + b= -2 (hence b = -1 ± 1^5^, 
2Re(t) = 1 and (from (^) 2Re{tb) = -35 - 5 gives 2t 
obtain — to? + bta — 56 = 0. The other two roots of this polynomial are d and c 
(which to call d follows from 2 Re (6 + c — oh) = 5), and its Galois group over Q[5, b] 
is (3^) < Z^ 

Computation of A{Zq). We get as before 
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df. The generator 2 G Zg sends 



Moreover, (14. 5 p yields gb = cf and ac = be 
a e d ^ a ^ e ^ d a, b f c h^b f c ^ b, and g ^'g ^ g. 

Up to the relations gb = cf, ac = be = df, there are precisely 17 Zg orbits in 
b,b, c,c, f, f of degree < 3. Using (14. 2 p with g = 3 and (14. 3 p when 3 divides both g 
and h, we can compute these 17 orbit sums Tr as linear expressions in s := Trfe = 
2Re(6+c+/). For example, Tr ^ = 1-s, Tr (bg) = 3(s-5), Tr (6^) = 25+l + s(2-5). 
But TrftTr^ = Tr (bg) + Tr (bg), so + (4 - 5)s + 1 - 5 = 0. Choose one of these 2 
solutions for s (they both yield A G ^(Zg)). 

From g + g = 1 — swe get g = {1 — s ± a/(5 — 6)s — 35)/2. Again select one 
of these g (they are complex conjugates). Let t := b + c + f. Then t + t = s and 
tg + tg = 3(s — 5), which fixes t = s/2 ± (35 — 6 — 5s — s) a/(5 — 6)s — 35/36 (use the 
same sign as in g). Then bc + bf + cf = tg and 6c/ = 51), so 6, c, / are the 3 roots of 
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— tx^ + gtx — 6g = 0. Let b be any of these roots; which to call c and / is fixed by 
the relation Tr (b'^c) = QS + 3 + {8 - 6)s. 

Now (H^ at (^f, h) = (1, 1) and (4,2) gives + va/u + 6u/u = for {u,v) = 
(b + (P/ f — 1,1 + be + be) and (/ + c — be/ 6,6 — f — f) respectively; subtracting these 
expresses a rationally in terms of b, e, f. Then d = af /e and e = be/ a. These values 
oi a, . . . ,g will solve (I4.1l) - (l4.3p . There are 2 choices for s, 2 for g, and 3 for b. The 
Zg symmetry accounts for this g and b ambiguity, leaving 2 inequivalent solutions. 

Computation of ^(Z|). Write the (l,0)-row of the matrix {6 — 1)74(^3) as 
{—l,—l,t,u,v,w,x,y,z). (14. 5 p with g = h = (1,0) and k = (0,1), (1,1) gives 
yu = xv = zw =: 6a for some complex number a with |a| = 1. From (14. 5p with 
g = (1,1), h = (1,0), A; = (1,2) yields 6t = xyz, and applying the Zg-automorphism 
(z, j) (-)■ (z, — j) to this we get 6t = uvw. Comparing these gives a^ = 1. 
Putting g = {0,l),h= (1, 0) into (g^D yields 

I - S = {au + au - l){av + av - 1) . (4.11) 

Hitting this with the Z|-automorphism ^ {i + j,j) gives the same equation 

with w replacing u. Thus au + au — 1 = av + av — 1 and (14. lip says there is a real 
number whose square is negative. This impossibility establishes part (f). 

Determining the modular data. Consider first K = and fix A G ^(Zs). Let 
Ui,Cgi^ be some solution to (14. 6p . Then we learn from the above proof that any 
automorphism a G Gal(Q/Q[5]) corresponds to an £ G Z3 in such a way that cAg^h = 
A.£g/h- Certainly a sends cui to some other root of unity of the same order d, and 
moreover auJi = uJi^ (since a commutes with complex conjugate in cyclotomic fields). 
Hence a sends the solution ui,Cgi^ of (14. 6 p to another solution uit,C^Ji^ given by 
^{^g,h) — ^eg,ih- Incidentally, (14. 7p then implies the Galois actions 

cr('S'j,,,0;,) = 5*0,^ _Bj,^ , (7(Ta;^aJ = T^i^^i^ 

(the second equation is simply the statement aui = ujit; to see the first, compute 
instead aS). This is stronger than the usual Galois action (II. 4p in modular data, 
though of course it is compatible with P^Hg^. 

Since 6 G Q[^^], when is coprime to the order d of uji any automorphism in 
Gal(Q[^rf]/Q) will lift to a a fixing 6 (when gcd(rf,/i) > 1, only half will). So if 
(14. 6 p can be solved for a di\i root of unity wi, for d coprime to then it can be 
solved for all di\i roots of 1. That is, we would have at least (f){d) solutions. As 
long as d doesn't divide v, none of these are redundant. But Izumi tells us there are 
precisely m nonredundant (independent) solutions to (14. 6 p (corresponding to the m 
half-braidings, giving rise to the m primaries of type 0/). 

Hence any solution to (14. 6 P involves Ui of order d where (l){d) < m (for d coprime 
to n) or (f){d) < 2m (otherwise). There is thus a finite number of possibilities for ui. 
A little thought reduces the number further — e.g. when gcd{d, //) = 1 it suffices 
to consider only ui = ^d- We find that (for u = 3) only 12 possibilities for ui need 
be considered, namely ±1, ±^3, 1, ±^7, ±^9, ±^13, —^13- To rule out these possibilities. 
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we compute determinants. In particular, for ui 7^ 1,^3, we show that the hnear 
system (14. 6 p is inconsistent by showing the (z/^ + 1) x (z/^ + 1) augmented matrix 
(formed from the coefficients and constant vector) has nonzero determinant. To 
show U[ = 1 (respectively ui = ^3) has at most 3 (resp. 1) independent solutions, 
we deleted 3 (resp. 1) equations and 2 (resp. 0) variables from (14. 6 p so that the 
resulting reduced coefficient matrix has nonzero determinant. We evaluated these 
determinants numerically using Maple. If some determinant was fairly close to 0, we 
found some Galois associate of ui with large determinant (i.e. of order 1 or higher). 

In this way we show that Ui must be one of the values . We can now determine 
the modular data, rigourously identifying it with P°Hg^, as follows. Numerically 
compute the solution to (14. 6 p and plug it into (14. 7p to obtain an estimate S'. We find 
that S' matches S to within 2 x 10"^, well within the 0.015 required by (14. 8p . 

This argument for K = and Z7 is identical. For them exactly 14 and 22, 
respectively, values of ui need to be considered. The error |S" — 5*1 is less than 
2 X 10"^ and 5 x 10^® respectively, again well within (14. Sp . 

The modular data calculation for = Zg is more complicated for two reasons, 
both related to the compositeness oi fi = 5-17. Fix A G A{1ig) (the result for the other 
A is automatic by Galois considerations). First, an automorphism a will permute 
the entries of A iff a fixes both 6 and s = Tr6 (incidentally, Gal(Q[s]/Q) = so s 
lies in a cyclotomic extension of Q). Hence potentially only 1/4 of Gal(Q[^(i]/Q) lifts 
to such a. This means more possibilities for ui to be eliminated (218 possibilities is 
overkill), and we find that the T matrix for A exactly matches that of either V^B.gg 
or V'^'^^'^'^^Ugg. More significantly, because of equalities like coi = wig, there won't 
always be unique solutions to (14.61) . and (14. 7p is ambiguous. This means we cannot 
identify S without more work. QED 

The key difference between K = and Z3 is the size of their automorphism 
groups (6 vrs 48). Nonuniqueness for K = Zg arises because /i = 85 is composite: 
more generally, where k > 1 is the number of distinct prime divisors of fi, there 
will be at least (2) inequivalent A and indeed inequivalent subfactors, realising the 
inequivalent modular data of type P°''^Hg^ defined in Section 3.3. 

This nonuniqueness of subfactors with identical principal graphs is not a surprise: 
for a simple example, the principal graph fails to uniquely identify the subfactor even 
in index < 4 (two subfactors realise each of the graphs Eq, Eg although one is the 
opposite of the other) whilst in index 4 the affine graphs ^2n-5 orbifold Dn^ 

both have n — 2 inequivalent subfactors [39] . 

Exact expressions for the matrices A e A{K) appear in the proof of Theorem 5. 
Numerical estimates for these A may also be of interest. A convenient way to express 
this, for K cyclic, is in terms of j2, is, • • • ,jn+i G for 0<g<h<uwe have 
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where ji = and jn+i+i = jn+i + jn — jn-i for 1 < i < n (see Lemma 7.3 of [38]). 

jf^ ^1.292076; 
of ~ (0.1846862, 1.5984702) ; 

(jf^jf^jr) -(2.471228,0.51685555,0.2137724); 

(jf \ • • -JP) ~ (2.396976693,2.079251103, -0.2079168419, -2.508673987) ; 
(jf . . . , jf^') ^ (-2.364737070, 1.031057162, 1.569692175,0.3383837765) . 

j{9) corresponds to I'^Hgg and j^^^' to ©''''"^^'^''j'Hgg. 

It should be possible to identify A{Z^) for the next couple z/'s, using the method 
of = 9 described in the proof. To probe the answer, we used Maple to obtain 
numerical solutions to fl4.ip -( H73l) . We don't have any proof that these correspond 
to actual solutions (although the numerics are convincing), and we certainly have 
no proof that there are no other solutions for these i>. The inequivalent numerical 
solutions we obtained for K = Z^, 11 < < 19, are: 

,(11) 



,(iiy 

,(13) 
,(15) 
,(17) 
,(19) 



0.9996507, 2.7258434, -0.5714203, -1.7797340, 1.2675985) , 
-2.6444397, -1.7629598, -2.6444440, 2.7572657, 0.1128260) ; 
-3.1050384, 0.5993399, -0.111708, -0.969766, 1.336848, 1.00483129) ; 
-1.0777623, -.7748018, -2.171863, -1.6068402, -.257508, 2.092502, .72289565) ; 
-1.466074, .291489, 3.130735, -2.693185, 1.398153, -.611938, -1.667078, -1.7548 
-2.677465, 1.088972, -.899442, .015448, -1.240928, -.493394, 1.839879, 

- 1.525884, -2.084374) ; 

^ (.896858, -.882585, -2.369855, -1.873294, -1.711620, -.119360, 2.972018, 

- 2.460652, .041334) . 

For u = 13, 15, 17 respectively, /i = 173, 229, 293 is prime. In these cases, choosing 
ui = ^'^^ yields a unique solution to fl4.6p (according to Maple) and plugging into 
fl4.7p yields a very close approximation (of order 10~^ or so) S' to the S matrix of 
P°Hg^. With high confidence we can expect j'^'^^ for u = 13, 15, 17 to describe a 
subfactor of Izumi type Zi, reahsing the modular data P^Hg^,. 

For z/ = 11,19 respectively, /x = 5^,5-73 are composite and the situation is 
more subtle. We find (up to the numerical accuracy of Maple) that for these 
uji = in fl4.6p yields the correct multiplicity of solutions, and thus the correspond- 
ing T matrix agrees with that of V'^Hg^. Likewise, the solution j^^^^' corresponds to 
-po,{5,73}jjg^^ As with u = 9, the existence of higher multiplicities means we cannot 
determine 5* unambiguously from (14. 7p for these 3 As. We expect j^^^^' 
to correspond to actual solutions of (14. 10 - 04. 3p . hence subfactors of Izumi type 
Zii,Zi9,Zi9 respectively, with modular data P^Hg^^, "D^Hg^^g, P°''!^^''''^^Hg]^g resp. 
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The only unexpected solution here is j^^^^'. Following the method described in 
the proof of Theorem 5, we obtain the corresponding T matrix: for each 1 < / < m, 

- \ ^5mi^ otherwise ' ^ ^ 

As before, S is not uniquely determined by (14.71) . We do not have a guess for what 
this S matrix is, though its entries Sd^^j,i, should lie in Q[cos(27r/25)]. Just as extra 
solutions like occur whenever distinct primes divide /i, we would expect 

that solutions like j^^^^' occur whenever a nontrivial prime power divides /i. 

On the basis of these observations, it is tempting to make the following conjecture. 

Conjecture 1. The modular data P^'^Hg^^ for any v = 2n+l and any set P (possibly 
empty) consisting of an even number of distinct prime divisors of fi, is realised by 
the even part of the quantum double of a subf actor of Izumi type Zy. All A G •A{K) 
should correspond to modular data obeying T'^t^'^ = /. 

We expect T^l^^^ = 1 for reasons clear from the proof of Theorem 5. Existence is 
accumulating suggesting that (up to equivalence) there is a unique matrix A G A{Z^), 
and a unique subfactor of Izumi type Zj,, iff = z/^+4 is prime. It would be interesting 
to see if A{K) can be nonempty when K is not cyclic. 

Question 5. Interpret the Zj, twist in P'^Hg^ in a von Neumann algebra formalism. 

We would expect this to be in terms of twisted systems of endomorphisms, or 
twisted fusion categories, as what happens for twisted doubles of finite groups. 



4.2 Principal graphs, a-induction etc for Izumi's subfactors 

Consider a subfactor C M of Izumi type K, an abelian group of odd order i/, 
defined by some matrix A G A{K), with A^-A^ and M-M systems A and A respec- 
tively. We argued last subsection that the modular data of the double T>A and VA 
is often (always?) the untwisted Haagerup-Izumi modular data V^Hg^^ or some gen- 
eralisation thereof. This subsection describes the data associated to this subfactor, 
generalising Section 2.2. For the special case of the Haagerup subfactor, i.e. K = Z3, 
much of this was computed in [5], though without appreciating the underlying Z3 
structure which organises everything. 

Write K for the inclusion A^ C M. Write K' for a set of representatives of the 
equivalence classes {K\0)/±; then g & K' labels the conjugacy classes \g\ := {g, —g} 
in the dihedral group _ft'xiZ2 (where Z2 acts by inverse). Recall 5 = {u + y/JI)/2, 
n = (z/ - l)/2, m = (z/2 + 3)/2, /i = z/^ + 4 = 2m + 1 and write A = y/1 + 6. 
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Figure 4. Principal graphs for Izumi's Z5 subfactor 

Theorem 6. Let N G M be a subfactor of Izumi type K as above. It has index 
5+1. The N-N system is A = {id, m, . . . , p, pw, . . . , pM,^_i}, and its sectors 
have fusions generated by 

WaWuh] = K+h] , KM = [p]Kg] , [p]' = [id] + J2[pug] . (4.13) 

The M-N system consists of endomorphisms nug (g G K) and 1% . The principal 
graph consists of v segments, each of length 3, sharing a common central vertex. The 
product M-N x N-N -+ M-N is 

[KUg][pUh] = W] + [KU_g+h] , [K'][pUg] = ^[kU h] + ( " , 

[nug][uh\ = [nug+h] , W][ug] = [k'] ■ (4.14) 
The product N-M x M-N N-N is: 

[KlQ[KUh] = [Uh-g] + [pUh+g] , [l^'][fiUg] = ^[pUh] , 

mK'] = J2K] + {u-l)^[pUg]. (4.15) 

The primaries of the double PA are 0, b, Oj, Cij, di for i G K' , jEK,l<l<m. 
The canonical endomorphism is 6'x)a = 1 + b + 2 with modular invariant = 

Ic/iq + chb + 2^^c/i|j,p. The alpha-inductions are a'^ = ax ® 1, a~ = 1 ® (q;x)°pp 
where, for all i G K' E K,l < I < m, 

[aa,] = 2[id] + J2^puj] , [ab] = [id] + ^^[pn,] , 
[acij] = [ui] + [m_,] + ^[pUj>] , [adi] = ^[pwj] . (4.16) 

Proof. In Lemma 7.1 and Theorem 7.2 of [38], we learn that A contains Ug and pUg for 
all g E K, and the fusions \Kk] = [id] + [p], [itc,][M/i] = [ug+h] and [p]^ = [id] + ^[pwg]- 



41 



By definition A consists of all irreducible sectors in any (kk)^, and from this we 
quickly obtain that A = {ug,pUg}. Because {up, Kp) = {Knp, p) = {{p + l)p,p) = 2 
and {up, nUg) = ((p + l)p,Ug) = Sg^o, we see that [up] — [k] := [k'] is an irreducible 
M-N sector distinct from any [nUg]. We obtain the fusion [n'Ug] = [k'] from the 
calculation {n'ug.K') = (^(p — l)ug,K{p — 1)) = ((p^ — l)%,p — 1) = 1. has 
dimension A^(5 — 1) = so will equal the sum of u sectors [pUg], perhaps with 
multiplicities, but i^-invariance now forces = '^[pug]. As above we see that 

the only sectors in nilltiY are [pwg], so we have exhausted all M-N sectors, and 
obtain the given principal graph. All remaining products are now easy to compute. 
For example, k'h' has dimension \^{5 — 1^ = u + {u — l)i'6, so by K-invariance this 

must be XIK] + (^^ - 1) 

At the end of Section 4 of [37] is explained how to recover the induction-restriction 
(i.e. dual principal) graph between A-B and A-A sectors, from the structure of 
Tube(A). In particular, we obtain the dual principal graph directly from the half- 
braidings as described in Proposition 8.2 of [38] — for instance 6 is read off from 
the half-braidings containing the identity. From the graph we read off the given 
alpha- induct ions. QED 

From this the dual principal graph for the double VA is the obvious analogue of 
the top graph of Figure 3. Note the resemblance between the fusions for A and the 
group algebra CD^ of = -ft'xiZ2, where Z2 acts on K by inverse. We also see that 
the sectors [ug], [pUg], [nug], and [k'] have statistical dimensions 1, 5, A and (5 — 1)A 
respectively. It is far from obvious (though of course an immediate consequence of 
Theorem 6) that the map a defined there is a ring homomorphism. 

Theorem 7. Assume the hypotheses and notation of Theorem 6. In addition, assume 
the fusions of the double system DA = {id, b, a^, Cjj, 0/} agree with those ofV^Hg^, 
and that we have the canonical endomorphism 9^-^ = l + b + ^aj + ^ Ch,o- Then 
we can identify the M-M system A with {id, a, bi,Ci}i^K' where Ci = KUjK = KU-iK. 
Its sectors obey the fusions 



where we set [cq] = = [bo] and write [c_j] = [q], [6„j] = [bi]. The dual principal 
graph is as in Figure J^: e.g. there is a (odd) central vertex, to which is attached [a], 
the [b,i\ 's and the [q] 's. The product M-M x M-N M-N is 



[ar = [id] + [a] + + , [a][k] = [a] + J2[b,] - [k] + ^^[c,] , 

H[c.] = [a] + J2\-^k] + [q] + 5^[cfc] , Mcj] = [a] + J^i^k] + ^^[cfc] , 
Mb,] = 6,,M + (1 - 6,,)[a] + - - + J2^c,] , 

Mcj] = 6i,,[id] + (1 + 6,,,)[a] + Y,h] + + + ^[Cfe] , 



[a][KUg] = [nug] + [k'] , = "^[nug] + (z/ - 1)[k'] , [b,i][KUg] = [k'] , 

[bi][K'] = {v - 2)[k!] + , = [ug+j\ + + [k'] , 




(4.17) 
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while the product M-N x N-M M-M is 



[KUg][KUh] = ([id] + [a])6g^h + [cg-h]{l - Sg,h) , [i^Ug][i^ = [a] + ^[k] + ^[q] , 

= [id] + (z. - l)[a] + (z. - 2) + • (4.18) 

The dual canonical modular invariant is Z = \ cHq + c/if, + X] '^^Oi + X] c/ic^ o P ■ Wfiiing 
a^ = ax®l and a~ = 1® i/ie alpha-inductions are 

[Sa,] = [id] + 2[a] + 5^ [6;,] - [h] + ^^[c,] , [ab] = [id] + [a] + ^^[6,] + ^[c,,] , 

[Sc,,,] = [id]5,,o + ([a] - [6,])(1 - 5,-o) + ^^[6^] + [q] + ^^[c,] , 

[ac),] = [a] + E[6.] + E[Q]. (4.19) 

Proof. From 6' we obtain immediately Z. The cardinality of A is determined from 
(7)7) = (^j ^) = XI ■^iQ = u + 1. Fix any / and define [id] = [So] (the identity sector 
of A), A = [a{di - E Ch„o + E B, = A + v [a{b - o^)], Ch = A + v [a{ch,o -h)]. 
We compute (using (ax, ay) = {9, xy) and the fusions of the double computed in 
Section 3.2) that (a(ai - b),a(ai/ - b)) = 1 + 5,,,/, (c^ - b, C/^/ - b) = 1 + 5h,h', 
(ttj - b, Cft - b) = -1, (t);, ai - b) = (3i, - b) = 0, (t);, 0;) = V. Hence the A, Bi, 
Ch are mutually orthogonal, each with norm z/, so they form a basis (over Q) for the 
vector space spanned by the sectors of A. 

Let [efc]o<fc<i/ be the irreducible sectors of A; choose [eo] = [id]. Then [A\, [Bi], [Ci\ 
are linear combinations (over Z) of these [cfc]: define A, L«, . . . , L^^") g Z'^+i by 

[A] = Y.^k [et] , m = [^] + ^ E [efc] , [C.] = [^] + E U'^'"'^ [ek] ■ Then i ■ 1 = 
u\ A-n^^ = -z/, ■ L« = 1 + 5k,i. From ([id],a„, - Sj,) = ([id],a,,_„ - St,) = 

([id],a5,) = we have Aq = L^^'^ = 0. Thus (reordering the [e^] if necessary) A = 
(0, si (z/ + t), . . . , s,y {u + t), st) and D'^'^ has components Lf'^ = sSi^y + SkSi^k, where 
s, Sk G {±1} and t G Z. This integer t satisfies u"^ = + u {u + t)"^, so t = —v. 

This means that [a] := ■^[^], and hence [fej] := ■^[Bi] and [cj] := ^[C*j], all being 
norm 1, are in A, possibly up to signs. That the signs are indeed all H-l's follows by 
writing say [SOi] in terms of [a], [6j], [cj]: the coefficients should be nonnegative. 

The fusions in A now reduce to the fusions in the double. The dual principal 
graph and the two remaining products follow by arguments as in Theorem 6. QED 

This guess for 9 is the natural one, matching what happens for the Haagerup 
subfactor (recall Section 2.2) and what we know about monomial modular invariants 
in Proposition 4, which appear to be severely constrained. The only example we know 
violating the hypothesis on the fusions of PA is the solution j'^^^^' of Section 4.1. 
The noncommutativity of A and commutativity of A arise because of the presence 
respectively absence of higher multiplicities in 9 resp. 9. The generalisation of Figure 
3 can be read off from the alpha-inductions listed in the Theorem. These imply the 
statistical dimensions dim(a) = 5, dim(6i) = 5 — 1 and dim(cj) = 5 + 1. 
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The dihedral group i^xiZ2 has one nontrivial 1-dimensional irrep and n 2- 
dimensional irreps. But for the Haagerup considered in Section 2.2, a corresponded to 
the 2- dimensional and b the 1-dimensional irreps, so K>iZ2 and A here are unrelated 
for u > 3, even after projecting away the Cj's. 

5 VOAs for Haagerup-Izumi modular data 
5.1 The Haagerup-dihedral diamond 

In Section 2.4 we saw "DHg as a mutation of VS^ and contained in a conformal 
subalgebra of T>1j3; now this becomes I^'^Hg^ is a mutation of V'^'^D^ and is contained 
in a conformal subalgebra of T>^^^. The branching rules for this embedding can be 
read off from the modular invariant fl3.1ip . In some ways (e.g. modular data) the 
Haagerup-Izumi closely resembles its sibling the dihedral group, and in other ways 
(e.g. VOA realisations) it more closely relates to its parent the cyclic group. 

The trivial case P*^Hg]^ is already interesting. In c = 8, its VOA can be realised 
as the tensor of the affine algebra VOA V(G'2,i) (which has c = 14/5), with the affine 
algebra VOA V(-F4 i) (which has c = 26/5). (Hence there is a VOA realisation for 
c any positive multiple of 8, by tensoring by copies of V{E^}j Call this V(Hg'j'). 
The VOA corresponding to PZi in c = 8 of course must be V(-E'8) = V{E^;y), and 
the containment corresponds to the conformal embedding G2, 1-^4,1 C -Es,!- The VOA 
corresponding to V-Di is the lattice VOA V{^AxE-j\ Although the latter is an orbifold 
of V{E^, V(G2, 1-^4,1) cannot be, since its modular data is not that of the (possibly 
twisted) quantum double of a finite group. This is an important lesson: we cannot 
expect the V^Hgj, VOAs to be orbifolds, by some subgroup of automorphisms, of 
V^Z^ VOAs. Hence Question 3. 

Another important observation from this trivial case = 1 is that the intersection 
of the holomorphic orbifold VDi with the u = \ Haagerup-Izumi VOA V(Hg'j') is 
itself a rational VOA. In particular, there are conformal embeddings V(i?3,iG'2.i) C 
V(E7,i) and V(Ai,i53,i) C V(F4,i) |56], and so both VD^ = V{AiEy) and V(Hg°) = 
V(G2, 1-^4,1) contain copies of the rational VOA V(Ai 16*2,1-83,1). Thus in this baby 
example at least, the dihedral and Haagerup-Izumi VOAs are commensurable in this 
sense. It is tempting to guess though that something similar happens for higher 
V. This is what we mean by the Haagerup-dihedral diamond: at the top is a VOA 
realising D'^Zy^ which contains as conformal subalgebras D'^'^D^ and X^'^Hg, and 
these contain a common conformal subalgebra (the analogue of V(Ai 16*2,153.1)). 

Conjecture 2. There is a rational VOA realisation of each V^Hg,^ for any suffi- 
ciently large central charge c (c a multiple of 8). 

As V increases, we will not be able to realise all P'^Hgj, — or even the D^Zy — as 
conformal subalgebras of V{E%). In other words, c > 8 will sometimes be necessary. 
This conjecture is open even for z/ = 3. To make it more plausible, and also to aid in 
the construction of the VOAs, in the following sections we supply all the information 
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needed to identify the possible character vectors. We saw in Section 2.4 that the 
Haagerup VOA would see the holomorphic orbifold by Z3 more directly than that by 
5*3, so for this reason we also supply the relevant data for V^Z^. 

We should repeat here the suggestion of Section 2.5 that a second approach to 
constructing VOAs realising V^Rg^ is to use the affine algebra V{Bm,2) in the coset 
construction. For reasons of space, we haven't explored this in this paper. 

5.2 Character vectors for the untwisted Haagerup double 

In this section we determine the matrices A and x for the untwisted Haagerup modular 
data P°Hg3, at all allowed values of central charge c (i.e. when 8|c). As explained 
in Section 1.4, from this all possible character vectors can be obtained. We illustrate 
this with c = 16 and c = 24 (for c = 8 see Section 2.4). 

5.2.1 Untwisted Haagerup at central charge c =24 16 

We found A, x for V'^B.g^ at c =24 8 in Section 2.4. Changing the central charge by 
a multiple of 8 leaves 5* unchanged but multiplies T by some 3rd root a; of 1. At the 
end of Section 1.4 we explain how this affects A, x- Taking the primaries in the usual 
order 0, b; a; Cq, Ci, C2; ^i, . . . , Og, we obtain for c =24 16 

A = diag(l/3, 1/3; 1/3; 1/3, 2/3, 0; 31/39, 7/39, 19/39, 28/39, -5/39, -2/39) , 





17 


155 


162 


162 


27 


729 


13 


286 


65 


13 


1001 


728 


155 


17 


162 


162 


27 


729 


-13 


-286 


-65 


-13 


-1001 


-728 




162 


162 


334 


-162 


-27 


-729 






















162 


162 


-162 


334 


-27 


-729 






















1215 


1215 


-1215 


-1215 


-76 


17496 






















9 


9 


-9 


-9 


6 


-12 






















1925 


-1925 














-21 


-6250 


-125 


50 


51331 


43175 




45 


-45 














-9 


58 


43 


12 


-45 


-297 




374 


-374 

















1650 


106 


-51 


-4250 


3927 




1288 


-1288 














23 


5313 


-759 


4 


39721 


-15432 


I 


1 


-1 














1 





-1 


1 


3 


1 


3 


-3 














3 


-9 


4 


-2 


6 


4 



From fll.24p the matrix S(r) is recursively computed. Any weakly holomorphic 
vector-valued modular function equals SP(J) for some polynomial vector P. 

Specialise now to a VOA character vector at c = 16. The corresponding poly- 
nomials Pi{x) are heavily constrained. For one thing, the vacuum character chQij) 
starts like lg~^/'^ + ■ ■ ■ and all other characters chi{T) start like a^g'^'"'^/^^ where 
hi > 0. This forces P = (1, 0, 0, 0, e, 0, /i, 0, j, 0, 0)"^ for constants e,g,h,j. So 
we've uniquely determined the hypothetical c = 16 Haagerup character vector, once 
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4 numbers (namely e,g, h,j) are specified. Its first few terms are 



/ choir) 

chair) = chcgir) 
c/lci (t) 
Chc2 (t) 
c/ioi(t) 
cho^ir) 
chjj^ir) 
chs^ir) 
\ ch^^^ir) I 
\ ch^eir) / 



+ (17 + 13g + 27e + 65h + 13j)q + (2013 + 594e + 4173/i + 260j + 91g)g2 H )\ 

' gi/^CllSS - 13g + 27e - 65h - 13i) + (21245 + 594e - 4173/1 - 260i - 91g)(?H ) > 

gl/3 ((-27e + 162) + (23247 - 594e)g + (705024 - 5967e)g2 ^ ) 

g-l/3 (e + (-76e + 1215)(? + (-1384e + 79704)q2 + (-11580e + 1886166)g3 H ) 

(6e + 9) + (5832 + 486e)g + (247131 + 5832e)g2 -\ 

^-8/39 + (-21g + 1925 - 125/1 + 50i)g + (-239g + 103631 - 4376/i + 506j)q^ H ) 

g7/39 ((43;j + 45 + 12j - 9g) + (4792/i + 258i + 10400 - 71g)g H ) 

g-20/39 (/i + (106/1 + 374 - 51j)g + (4980/1 - 634j + 34606 - 43g)g2 h ) 

g-il/39(j _,_ ^ 1288 + 23g - 759h)q + (-84j + 141c/ + 79583 - 25135h)q^ H 

g-5/39 + 1 + J + + (_856/i + 130i + 2704 + 31g)<? + ■ • ■ ) 
\ g-2/39 ((4/j 3 _ 2j + 3g) + (1582/1 - 95i + 3898 + 70g)q H ) 



172 


324 


324 


54 


1458 


162 


334 


-162 


-27 


-729 


162 


-162 


334 


-27 


-729 


1215 


-1215 


-1215 


-76 


17496 


9 


-9 


-9 


6 


-12 



Since the Fourier coefficients of a VOA character vector are nonnegative integers, 
we know each e,g, h,j G Z>o. In addition the g^/^-coefficient of cha forces e < 6 and 
the g^/^-coefficient of chi, then forces g + j + 5h < 24, so there are only finitely many 
possibilities. Other coefficients yield further bounds, e.g. h < 2, j < 11, g < 16. 

There are several other constraints. One is that VOA characters are linear 
combinations over Z>o of the Virasoro characters at that c for various values of 
conformal weight h (because the VOA carries a Virasoro representation). All 
these characters are well-known: c/i(i6,o) = q~'^^^{l + + + + ■■■) and 
ch{2A,h) = q'^^^^^il + q + 2q^ + 3q^ + 5q^ + ■ ■ ■) for h > 0. This isn't so useful 
here, but it is effective for the c = 24 case considered next. 

More important, Dong-Mason [18] prove that the conformal weight 1 part of a 
rational VOA V is a reductive Lie algebra (BiQi of central charge '^■kidimQi/{ki + 
h^) < c and dimension equal to the coefficient of g^"^/^^ c/^q. This gives a finite set 
of possible values for that coefficient, and knowing this reductive Lie algebra helps 
in constructing V. Each V-module is also a module for the corresponding affine Lie 
algebra, so the V-characters are linear combinations of the affine characters. 

This is as far as we'll carry the analysis here. 

The matrices A, x for with primaries in the order (00), (01) = (02), (10) = 

(20), (11) = (22), (12) = (21) are A = diag(l/3, 1/3, 1/3, 2/3, 0) and 
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At c = 16 there is only one free parameter, d: the possible character vectors are 

c/i(o,o) \ / 9-2/^(1 + (172 + 54d)g + (23258 + 1188d)q2 H ) 

ch(o,i) = c/i(o,2) = c/i{i,o) = c/i{2,o) I _ I (J^''^((162 - 27d) + (23247 - 594d)q + (705024 - 5967d)q^ H ) 

c/i(i_i) = c/i(2,2) I ~ I (j-i/3((i+ (1215 - 76(i)ij + (79704- 1384d)q2 H ) 

c/i(i,2) = c/i{2,i) / y (9 + Od) + (5832 + 486c/)g + (247131 + 5832(i)g2 H 

This implies d G Z>o and d < 6. This d equals the parameter e of V^Hg. 
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5.2.2 Untwisted Haagerup at central charge c =24 



Next turn to the Haagerup at 24|c. Finding its A, x is now routine: 

A = diag(0, 0, 1, 1, 1/3, 2/3, 6/13, 11/13, 2/13, 5/13, 7/13, 8/13) , 








1/2 


1/2 


18 


6 


6 


1 


10 


5 


4 


4 




-12 


1/2 


1/2 


18 


6 


-6 


-1 


-10 


-5 


-4 


-4 


65610 


65610 








-5832 


-243 




















65610 


65610 








-5832 


-243 




















729 


729 








-152 


54 




















8748 


8748 








2430 


-76 




















1716 


-1716 














-252 


-5 


-176 


176 


134 


100 


22451 


-22451 














-980 


8 


16556 


2464 


-96 


-385 


104 


-104 

















5 


44 


-56 


-15 


20 


910 


-910 














77 


8 


-847 


-32 


120 


-34 


3003 


-3003 














330 





-1540 


704 


-216 


55 


\ 5200 


-5200 














616 


-14 


3388 


-595 


132 


-20 



Its possible character vectors are found as before to be 

g-i + (-12 + 10i + 18e+ +6/ + 4fc + 5j + /i + 4« + 6fc) H \ 

(-lOi + 18e + £±^^ + 6/ - 4fc - 5i - h - 4Z - 6g) + • • • 

c + (65610 - 5832e - 243/)(j + (7164612 - 247131e - 2916/)g2 -\ 

d + (65610 - 5832e - 243/)g + (7164612 - 247131e - 2916/)q2 + . . . 

eg-2/3 _,_ (729 _ I52e + 54/)qi''3 ^ (370332 - 23236e + 1188/)g*/^ H 

/q-l/3 + (8748 - 76/ + 2430e)q2/3 + (1743039 - 1384/ + 159408e)(?5/3 + ■ ■ ■ 

gq-7/13 + (1716 _ 252g - 176i + 134fc + 100/ - 5h + nQj)q^/^^ H 

hq-^/^^ + (22451 + 8fc + 2464j + 16556i - 385/ - 96fc - 980g)g"/i3 H 

jg-ii/13 _|_ (;^Q4 _,_ „ i^f, + 5h- 56j + 20Z)g2/i3 ^ 

jg-8/13 _,_ (gj^o - 32j - 34/ - 847j + 77g + 8h + 120fc)g5/i3 _| 

fcq-6/13 _,_ (3003 + 704j - 1540j + 55/ - 216/c + 330g)q'^^^'-^ H 

lq-5/rA _|_ (5200 - 20/ + 3388i + 616g - 14h + 132/c - 595j)(?*/i3 H / 

for certain nonnegative integers c,d,e, f,g,h,i,j,k,l. By the usual arguments we 
find a finite number of possibilities. Likewise the matrices for V^Z^ are A = 
diag(0, 1,1,1/3,2/3) and 

(-12 11 36 ^2 \ 

65610 -5832 -243 \ 

65610 -5832 -243 

729 -152 54 I 

8748 2430 -76 / 

and we find a character vector depending on 4 bounded nonnegative integers. 

The most intriguing possibility is if these Haagerup and Z3 VOAs are both sub- 
VOAs of the Moonshine module. It is tempting to guess that this is the most natural 
home of the Haagerup. Then the character J of the Moonshine module will equal 
one of cho + cht, + 2cha, cho + ch^ + 2ch^, or ch^ + ch^, + cha + c/i^q, and we find that 
(c,d, /) is either (12,0,0), (0,12,0) or (0,0,1), and that {g,h,j,k,l) is either (1,2,0,0,0), 
(0,1,1,0,0), (0,2,0,1,0) or (0,1,0,0,1), and all other parameters vanish. 

5.3 Character vectors of the twisted Haagerup double 

The modular group representation for the w 7^ twist of the Haagerup decomposes 
into a sum of 3 irreps: the trivial one, a 4-dimensional one with kernel r(9), and 



/ cho 
cht, 
cha 
chcg 
c/ici 
chc2 
chj,-^ 

Chx)2 

chj)^ 
choi 
chj,^ 
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a 7-dimensional one with kernel r(13). Only the 4-dimensional irrep is new: it is 
handled by e.g. the root lattice A^. Actually, 6 different 4-dimensional irreps arise 
here, varying with the twist and central charge mod 24, but they are all obtained from 
the one in by some combination of 3rd roots of unity and taking the contragredient. 
The effects on A, x of both of these is discussed at the end of Section 1.4. 



5.3.1 The 1-twisted Haagerup at central charge c =24 8 

The matrices A,x for V^Z^ and V^Rg^ are respectively diag(2/3, 2/3, 7/9, 1/9, 4/9), 





80 168 


54 


9504 


1078 \ 




















84 164 


-27 


-4752 


-539 


















( 


126 -126 


8 


-17248 


1375 




















9 -9 


-8 


96 




10 1 




















36 -36 


20 


392 




-340 / 


















dia 


g(-l/3, 2/3, 2/3, 7/9,1 


/9, 4/9, 5/39, 20/39, 32/39, 41/39, 8/39, 11/39) , 






/ 


-42/5 


1 














-17/5 


-13/5 


3/5 


1/5 


-7 


14/5 


\ 


209456/5 


80 


84 


81 


14256 


1620 


-19054/5 


-2366/5 


166/5 


77/5 


-2814 


-5812/5 




228912/5 


168 


164 - 


-81 


-14256 -1620 


-1428/5 


-1092/5 


252/5 


84/5 


-588 


1176/5 






176556/5 


84 


-42 


8 


-17248 1375 


-714/5 


-546/5 


126/5 


42/5 


-294 


588/5 






114/5 


6 


-3 


-8 


96 


10 


-51/5 


-39/5 


9/5 


3/5 


-21 


42/5 






11391/5 


24 


-12 


20 


392 


-340 


-204/5 


-156/5 


36/5 


12/5 


-84 


168/5 






546/5 

















-559/5 


-161/5 


16/5 


7/5 


14 


418/5 






52611/5 

















-8134/5 


-86/5 


206/5 


42/5 


-486 


-4617/5 






600457/5 

















-12733/5 


3388/5 


412/5 


119/5 


-6868 


24871/5 






2851862/5 

















100947/5 


6118/5 


437/5 


229/5 


31027 


-49514/5 




\ 


1729/5 

















649/5 


-174/5 


14/5 


13/5 


-238 


437/5 


/ 


5252/5 

















2002/5 


-567/5 


77/5 


14/5 


27 


76/5 



At c = 8, T>^1^3 is realised by the lattice VOA V{As). "D+'^^s is realised by the affine 
algebra VOA V(i?4,2)- Both of these are conformal embeddings. 



5.3.2 The 1-twisted Haagerup at central charge c =24 16 

The matrices A, x for V^Z^ and V^B.g^ are respectively diag(l/3, 1/3, 4/9, 7/9, 1/9), 



166 


330 


198 


18 


900 


165 


331 


-99 


-9 


-450 


351 


-351 


56 


-25 


1625 


2079 


-2079 


-1694 


53 


3146 


27 


-27 


28 


1 


-102 



diag(-2/3, 1/3, 1/3, 4/9, 7/9, 1/9, 31/39, 7/39, 19/39, 28/39, 34/39, -2/39) , 
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-1 
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-1 
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-1 


6732 


166 


165 


297 


27 


1350 


-165 


-286 


87 


-165 


152 


-880 


7359 


330 


331 


-297 


-27 


-1350 


-165 





165 


-165 


165 


-165 


8199 


234 


-117 


56 


-25 


1625 


-117 





117 


-117 


117 


-117 


173727 


1386 


-693 


-1694 


53 


3146 


-693 





693 


-693 


693 


-693 


126 


18 


-9 


28 


1 


-102 


-9 





9 


-9 


9 


-9 


498225 

















-1946 


-6250 


1800 


-1875 


1925 


41250 


858 

















-54 


58 


88 


-33 


45 


-342 


31603 

















-374 


1650 


480 


-425 


374 


3553 


263120 

















-1265 


5313 


529 


-1284 


1288 


-16720 


920556 

















-2673 


-2025 


1848 


-2574 


2704 


22572 


13 




















-9 
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-5 
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1 
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5.3.3 The 1-twisted Haagerup at central charge c =24 



The matrices A, x for ^^^^3 and V^Yig^ are respectively diag(l, 0, 1/9, 4/9, 7/9), 






131274 


185328 


13770 


324 


1 


6 


-135 


-54 


-27 





-27 


24 


-35 


10 





-594 


-2002 


290 


11 





-5967 


14144 


340 


-64 



diag(0, 1, 0, 1/9, 4/9, 7/9, 6/13, 11/13, 2/13, 5/13, 7/13, 8/13) , 



/ -12 
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20 


2 


10 


8 


12 


8 


60033 





65637 


92664 


6885 


162 


-18997 


-16 


-3200 


-350 


-1617 


-792 


-18 
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-135 


-54 


-27 


20 
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10 
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12 


8 


27 





-27 


24 


-35 


10 
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-594 


-2002 


290 


11 




















5967 





-5967 


14144 


340 


-64 




















1716 

















-252 


-5 


-176 


176 


134 


100 


22451 

















-980 


8 


16556 


2464 


-96 


-385 


104 




















5 


44 


-56 


-15 


20 


910 

















77 


8 


-847 


-32 


120 


-34 


3003 

















330 





-1540 


704 


-216 


55 


V 5200 

















616 


-14 


3388 


-595 


132 


-20 



5.3.4 The 2-twisted Haagerup at central charge c =24 8 

The matrices A, x for ^^^^3 and V'^Yig^ are respectively diag(— 1/3, 2/3, 8/9, 2/9, 5/9), 
2 1 1/3 -1/3 -2/3 \ 

46683 248 26 -836 -133 \ 

706401 156 2584 475 , 

2187 14 172 -77 / 

56862 65 -2431 -74 / 



diag(-l/3, -1/3, 2/3, 8/9, 2/9, 5/9, 5/39, 20/39, 32/39, 41/39, 8/39, 11/39) , 



/ -16/5 
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1/6 


-1/6 


-1/3 


-17/10 


-13/10 


3/10 


1/10 


-7/2 


7/5 


26/5 


-16/5 
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1/6 


-1/6 


-1/3 


17/10 


13/10 


-3/10 


-1/10 


7/2 


-7/5 


46683 


46683 
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26 


-836 


-133 
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706401 





156 


2584 
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-2431 


-74 
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-546/5 














-559/5 


-161/5 


16/5 


7/5 


14 


418/5 


52611/5 


-52611/5 














-8134/5 


-86/5 


206/5 


42/5 


-486 


-4617/5 


600457/5 


-600457/5 














-12733/5 


3388/5 


412/5 


119/5 


-6868 


24871/5 


2851862/5 


-2851862/5 














100947/5 


6118/5 


437/5 


229/5 


31027 


-49514/5 


1729/5 


-1729/5 














649/5 


-174/5 


14/5 


13/5 


-238 


437/5 


\ 5252/5 


-5252/5 














2002/5 


-567/5 


77/5 


14/5 


27 


76/5 
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5.3.5 The 2-twisted Haagerup at central charge c =24 16 

The matrices A, x for P^Za and V'^Rg^ are respectively diag(l/3, 1/3,5/9,8/9,2/9), 



v 



160 


336 


32/3 


2/3 


196/3 
















168 


328 


-16/3 


-1/3 


-98/3 
















5832 


-5832 


-272 


2 


1925 
















32076 


-32076 


196 


12 


-15092 
















729 


-729 


40 


-4 


28 
















g(-2/3, 1/3,1 


/3, 5/9, 8/9, 2/9, 31/39, 7/39, 19/39, 2t 


V39, 34/39,- 


2/39) , 




3 


1 














-1 





1 


-1 


1 


-1 


6717 


160 


168 


16 


1 


98 


-162 


-286 


84 


-162 


149 


-877 


7374 


336 


328 


-16 


-1 


-98 


-168 





168 


-168 


168 


-168 


221859 


3888 


-1944 


-272 


2 


1925 


-1944 





1944 


-1944 


1944 


-1944 


3788856 


21384 


-10692 


196 


12 - 


-15092 


-10692 





10692 


-10692 


10692 


-10692 


5589 


486 


-243 


40 


-4 


28 


-243 





243 


-243 


243 


-243 


498225 

















-1946 


-6250 


1800 


-1875 


1925 


41250 


858 

















-54 


58 


88 


-33 


45 


-342 


31603 

















-374 


1650 


480 


-425 


374 


3553 


263120 

















-1265 


5313 


529 


-1284 


1288 


-16720 


920556 

















-2673 


-2025 


1848 


-2574 


2704 


22572 


13 




















-9 


7 


-5 


3 


1 



5.3.6 The 2-twisted Haagerup at central charge c =24 

The matrices A, x for V^Z^ and V^Rg^ are respectively diag(l, 0, 8/9, 5/9, 2/9), 



131274 


6 


528 


9282 


-3 


-1 


-7 


-8 


-104247 


3 


-1001 


12376 


-12393 


-7 


232 


476 


-729 


5 


22 


-224 



diag(0, 0, 1, 2/9, 5/9, 8/9, 6/13, 11/13, 2/13, 5/13, 7/13, 8/13) , 



-15/2 


9/2 


1/2 


4 


7/2 


1/2 


6 


1 


10 


5 


4 


4 


9/2 


-15/2 


1/2 


4 


7/2 


1/2 


-6 


-1 


-10 


-5 


-4 


-4 


65637 


65637 





-4641 


-264 


-3 




















729 


729 





-224 


22 


5 




















12393 


12393 





476 


232 


-7 




















104247 


104247 





12376 


-1001 


3 




















1716 


-1716 














-252 


-5 


-176 


176 


134 


100 


22451 


-22451 














-980 


8 


16556 


2464 


-96 


-385 


104 


-104 

















5 


44 


-56 


-15 


20 


910 


-910 














77 


8 


-847 


-32 


120 


-34 


3003 


-3003 














330 





-1540 


704 


-216 


55 


5200 


-5200 














616 


-14 


3388 


-595 


132 


-20 
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